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^! Abstract 
D 

' We present a study of a class of effective actions which, show typical axion-like interactions, and 

of their possible effects at the Large Hadron Collider. One important feature of these models is the 
^ . presence of one pseudoscalar which is a generalization of the Peccei-Quinn axion. This can be very 

CN , light and very weakly coupled, with a mass which is unrelated to its couplings to the gauge fields, 

described by Wess Zumino interactions. We discuss two independent realizations of these models, 
' one derived from the theory of intersecting branes and the second one obtained by decoupling one 

l/^ ' chiral fermion per generation (one right-handed neutrino) from an anomaly-free mother theory. The 

. key features of this second realization are illustrated using a simple example. Charge assignments 

, of intersecting branes can be easily reproduced by the chiral decoupling approach, which remains 

more general at the level of the solution of its anomaly equations. Using considerations based on 
its lifetime, we show that in brane models the axion can be dark matter only if its mass is ultralight 
^ ■ (- 10"'' eV), while in the case of fermion decoupling it can reach the GeV region, due to the 



absence of fermion couplings between the heavy Higgs and the light fermion spectrum. For a GeV 
axion derived from brane models we present a detailed discussion of its production rates at the 
LHC. 



1 



1 Introduction 



The study of possible signatures of string/brane theory at lower energy has achieved a significant 
strength with the development, in the last few years, of several extensions of the Standard Model 
(SM) formulated in scenarios with intersecting branes and large extra dimensions [B El El HI |6] , 
which are characterized by quite distinct features compared to other constructions, such as those based 
on more traditional anomaly- free supersymmetric formulations. The latter include specific theories 
like the MSSM but also its further variants such as its next-to-minimal (nMSSM,NMSSM) extensions, 
eventually with the inclusion of a gauge structure enlarged by an extra anomaly- free ^7(1) gauge 
symmetry (USSM) [7] (see [8j for an overview). 

On the other hand, since anomalous C/(l)'s are naturally produced in geometrical compactifications 
and are an important aspect of brane models, the search for possible signatures of string theory 
has necessarily to take into consideration the peculiarities of these anomalous extensions, which are 
characterized by anomalous extra neutral currents, contact interactions of Chern-Simons form at 
trilinear gauge level [SJ] [TU] [TT] and several axions of Stiickelberg type. Supersymmetric extensions of 
these classes of models have also been investigated recently [121 13] • 

One of the most demanding feature of these formulations, in regard to possible experimental 
searches, is to clarify the role of gauge anomalies on a substantial sector of collider phenomenology, 
from precision measurements of leptoproduction to double prompt-photon production [T3], just to 
mention a few processes. These and many more are all affected by the new anomalous trilinear gauge 
vertices |15) which appear in these models, although their studies are expected to be quite difficult 
experimentally. 

In fact, the limited accuracy of hadron colliders might reduce the expectations in regard to the 
possible experimental identification of subtle effects due to the mechanism(s) which underline the can- 
cellation of the gauge anomalies. Nevertheless, the presence of an axion-like particle in the spectra of 
these theories is an important feature of intersecting brane models, which represents a serious depar- 
ture from the typical anomaly-free formulations - both for supersymmetric and non-super symmetric 
models - and provides a natural justification for a light pseudoscalar state. 

The higher perturbative order at which these effects start to appear in the perturbative expansion 
and the limitations of the parton model description seem to indicate that the analysis of anomalous 
effects are more likely to be the goal of a linear collider rather than that of the LHC, nevertheless the 
signatures of new physics are manyfold and are not limited to collider physics, but have remarkable 
implications also in astroparticle physics and cosmology. 

Among the aspects that can be addressed within these new formulations are those related to the 
flavour sector and the connection between these constructions and the traditional solutions of the 
strong CP-problem, previously addressed with the help of global U{1) symmetries, such as in the 
invisible axion model [16\ \T7\ [T8\ \T9\ [20] . We recall that studies of the flavour sector of the SM in 
the presence of gauged anomalous C/(l)'s are not new, having been used in the past in a variety of 
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cases, for example in the construction of realistic scenarios for neutrino mixing [21]. At the same time, 
the study of axion-like particles is at the center of new important proposals for their detection which 
are now under an intense investigation at DESY 122\ I23j . Other interesting proposals consider the 
possible implications of axion-like particles in the propagation of gamma rays [24J. We believe that 
these motivations are sufficient to justify generalized searches of pseudoscalars as a possible solution 
of the dark matter problem. At the same time anomalous gauge interactions, in combination with 
quantum gravitational effects, show puzzling features, due to the presence of phantom fields [251 ES] 
in the local formulation of the trace anomaly [2^ which deserve a closer look. 

1.1 An axion with independent gauge couplings and mass 

An axion-like particle is characterized by the usual pseudoscalar couplings to the gauge fields (the 
bFF term, where b is an axion) but has a mass which is unrelated to its coupling. Different mass 
ranges for the axion have quite different implications at a phenomenological level. For instance, for a 
very light axion (~ 10~^ — 10~^ eV), as for the PQ case, the pseudoscalar can mix with the photon 
and can generate, in the presence of background galactic magnetic fields, the usual phenomena of 
birifringence and dichroism for light propagation ^28j, with important effects at astrophysical level 
\24:\ I29j and other experimental signatures |22[ I23| . The optical activity of the intergalactic medium 
due to the presence of background axions, also in this generalized case, is essentially caused by the 
bFF coupling in the equations of motion of the lagrangean [3D] [31] . 

In the invisible axion model astrophysical arguments bound the mass of the axion (and its inter- 
action to the gauge fields) requiring its suppression by a large scale /. All the axion couplings and 
the axion mass 

m PS 6 • lO^^eV-pGeV (1) 

are inversely proportional to /, where / is arbitrary (/ ~ 10^ GeV experimentally) and makes the 
axion, indeed, very light. In general, a very light axion, being a quasi-goldstone mode of a global 
symmetry, is produced copiously at the center of the sun and escapes after its production, with a 
mean free path which is larger than the radius of the sun. The failure by existing ground-based 
helioscopes to detect this particle in a detector of Sikivie type [32] [221 El] has been used to bound 
its mass and its interaction with the gauge fields. The bound can be evaded if the axion has a mass 
larger than the temperature at the center of the sun, since in this case would not be produced at 
its center, mass which is not allowed for the invisible axion according to current constraints. For a 
very light axion interesting effects are allowed, such as its non-relativistic decoupling, since its average 
momentum at the QCD phase transition is not of the order of the associated temperature, which is 
in the GeV range, but of the Hubble expansion rate (3 • 10~^ eV), and the formation of Bose-Einstein 
condensates [35] . 

As we have mentioned, the gauging of the axionic symmetries can lift the typical constraints of 
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the invisible axion model, allowing a wider parameter space, which is the main motivation for our 
study. In principle, in the extensions that we consider, this pseudoscalar can be very light, while 
its gauge interactions can be suppressed by a scale which is given by the mass of the lightest extra 
Z' present in the neutral sector of these models. For this reason, these types of pseudoscalars are 
naturally associated to the neutral current sector, with new implications at the level of the trilinear 
gauge interactions. 

So far, two models have been developed in which the structure of the effective action allows 
a physical axion: the MLSOM (the Minimal Low-Scale Orientifold Model) [36] and the USSM-A 
[371 113| . the first being a non-super symmetric model, the second a supersymmetric one. In the first 
model, motivated by a construction based on intersecting branes, the scalar sector involves beside the 
Stiickelberg axions, 2 Higgs doublets. At the same time, the gauge structure of the Standard Model 
is corrected by the presence of extra neutral currents due to the extra C/(l). 

To date, a detailed analysis of these models is contained in [38j, worked out for a single extra C/(l). 
In the supersymmetric case the presence of a physical axion is guaranteed if the superpotential allows 
extra superfields which are singlet respect to the Standard Model but are charged under the anomalous 
C/(l)'s. The field content of the superpotential of the nMSSM is sufficient to have a physical axion in 
the spectrum [371 [T3]. 

1.2 Gauged axions 

The gauging of axionic symmetries is realized in the low energy effective lagrangeans by introducing 
(shifting) axions, one for each anomalous C/(l) present in the gauge structure of a given model. These 
are accompanied by Wess-Zumino terms in order to restore the gauge invariance of the theory due 
to the chiral anomalies present in these constructions. These axions (Stiickelberg axions) are not all 
physical fields. In fact, the only physical axion, called the "axi-Higgs" in [36], is identified in the 
CP-odd sector of the scalars by a joint analysis of the potential and of the bilinear mixing terms 
(Bidbi) generated by the Stiickelberg mass terms which are present for each anomalous U{1). We 
are going to summarize below the scalar of the scalar potentials which allow a physical axion in the 
spectrum, either massless or massive. Since the mass of this particle is expected to obtain small 
non-perturbative corrections due to the instanton vacuum, as for the invisible axion, these small 
corrections are described by extra terms in the scalar potential which are allowed by the symmetry. 
These terms make the physical axion part of the scalar potential, but their size remains, in the class 
of theories that we analyze, essentially unspecified. In supersymmetric models they are expected 
to correspond to non-holomorphic corrections to the superpotential [37] which involve directly the 
axion/axino superfield. The size of these corrections depends on the way the fundamental symmetry 
is broken, and the appearance of the axion in the scalar potential just parameterizes our ignorance of 
the fundamental mechanism which is responsible for these corrections. For this reason, we focus our 
analysis on several mass windows for this particle, although the most relevant mass range for collider 
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studies is the GeV region. 

1.3 Organization of this work 

The analysis presented in this work concerns the phenomenology of the axi-Higgs in anomalous abelian 
models with a single anomalous extra U{1) and in the non-supersymmetric case. The construction, 
therefore, is the one typical of the MLSOM, formulated in the context of intersecting branes. A 
similar analysis can be performed in the supersymmetric case, although it is more complex and will 
be presented elsewhere. Our analysis, however, is not limited to models of intersecting branes, but to 
the entire class of effective actions which are characterized by axion-like interactions at low energy, 
independently from their high energy completion. Typical charge embeddings of brane constructions, 
as we are going to show, can nevertheless be obtained in our approach starting from an anomaly-free 
spectrum and decoupling some chiral fermions. Some differences between the two realization remain, 
at phenomenological level, since the corresponding axion, in the case of decoupled fermions, does not 
couple to the light fermions which are part of the low-energy spectrum. 

Our motivations for working within this more general framework has been motivated by scenarios 
where a heavy fermion, for instance a right-handed neutrino, decouples from the low energy spectrum 
leaving one Stiickelberg axion (the phase of a Higgs field) in the effective lagrangean. We will come 
to discuss these points in more detail in one of the sections below. The different completions of these 
lagrangeans start differing at the level of operators whose mass dimensions is larger than 5, the five 
dimensional ones being the Wess-Zumino terms. 

After reviewing briefiy these models in order to make our analysis self-contained, we illustrate 
how their anomalous content can be obtained by requiring that only some of the anomaly equations 
are satisfied, taking as a starting point an anomaly-free chiral spectrum and decoupling some chiral 
fermions. Typical brane models such as the Madrid model [3] are obtained for a particular choice of the 
free charges allowed by the decoupling of the heavy chiral fermions and are just particular solutions 
of the anomaly equations. We then move towards a phenomenological analysis of the axi-Higgs in 
the MLSOM, selecting the GeV mass range for the axion. This region is the most promising one for 
collider studies of this particle, although in this range, as we are going to show, it is not long-lived. A 
GeV axion can be long lived, but must have suppressed couplings to the fermions of the low-energy 
spectrum, and one way of getting this lagrangean is via the mechanism of decoupling of heavy fermions 
(and of the radial excitations of the associated Higgs field) from the low energy theory. We show, 
using a simple toy model, how this can occur. 

Production and decay rates for this particle are studied for all the mass windows in the MLSOM 
for typical LHC searches. We give in an appendix a summary of the scalar sector of the lagrangean 
and the determination of coefficients of the Wess Zumino terms. We have also included a section where 
we present a discussion and a comparison of the effective action of intersecting brane models versus 
the analogous one obtained by decoupling a chiral fermion, illustrating briefly the origin of the various 
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operators left in the low energy formulation, with the axion interpreted as the phase of a second Higgs 
sector, partially decoupled from the 2 Higgs doublets included in the electroweak sector. 

2 The model: overview of its general structure 

We analyze a class of models characterized by a gauge structure of the form SU (3) x SU{2) x U (l)y x 
U{1)b-, defined in [38j, where the U{1)b gauge symmetry is anomalous and the corresponding gauge 
boson {B) undergoes mixing with the rest of the gauge bosons of the Standard Model. Details can 
be found in [38^ [38l I39j : here we just summarise the main features of this construction for which we 
will define rather general charge assignments. As we have already stressed, the reason for keeping 
our analysis quite general is motivated by the observation that effective actions of intersecting brane 
models are not uniquely identified. Various completions can generate the same low energy signatures, 
at least up to operators of dimension 5, which, for anomalous gauge theories, are the Wess-Zumino 
terms. These points will be illustrated in a section below, where we will solve the basic equations that 
characterize the charge assignments of the anomalous model, under some assumptions on the fermion 
spectrum which are essential in order to make our analysis concrete. 

2.1 The structure of the effective action 

The effective action has the structure given by 



where Sq is the classical action which is given in an appendix. It contains the usual gauge degrees 
of freedom of the Standard Model plus the extra anomalous gauge boson B which is already massive, 
before electroweak symmetry breaking, via a Stiickelberg mass term. The scalar potential is the 
maximal one permitted by the symmetry and allows electroweak symmetry breaking. The structure 
of the Yukawa sector Syuk is very close to that of the Standard Model. In one of the sections below 
we identify the fundamental physical degrees of freedom of this sector after electroweak symmetry 
breaking, which, in our analysis, is based on the choice of the largest potential allowed by the symmetry. 
The model is a canonical gauge theory with dimension-4 operators plus dimension 5 counterterms of 
Wess-Zumino type. 

In Eq. ([2|) the anomalous contributions coming from the 1-loop triangle diagrams involving abelian 
and non-abelian gauge interactions are summarized by the expression 



S 



So + Svuk + San + SwZ + SqS 



(2) 



5, 



an 



1^{TbwwBWW) + ^{TbggBGG) + ^{TbbbBBB) 
+^{TbyyBYY) + ^{TybbYBB), 



(3) 
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where the symbols () denote integration. For instance, the anomalous contributions in configuration 
space are given explicitly by 

{TbwwBWW) = j dxdy dzTX^{z,x,y)B^{z)Wi'{x)W^ (y) (4) 

and so on, where Tbww denotes the anomalous triangle diagram with one B field and two W^s external 
gauge lines. The gluons are denoted by G. 

In the same notations the Wess Zumino (WZ) counterterms are given by 

Swz = ^{bFBAFB) + ^{bFYAFY) + ^{bFYAFB) 

+ ^(6rr[F^ A F^]) + ^{bTr[F^ A F^]), (5) 

while the gauge dependent CS abelian and non abelian counterterms [10] needed to cancel the mixed 
anomalies involving a B line with any other gauge interaction of the SM take the form 

Scs = +di{BY AFy) +d2{YB AFb) 

+c,{e^^'^P"B^C^,l!f)) + C2{e^^P'^B^Cl^}^)), (6) 



with the non-abelian CS forms given by 



f<SU{2) 



K [^Z^P + 1^2 e^'^'WlWl^ + cyclic 



Gl[Fl,p^-g^r'^GlG^^\ acyclic 



(7) 
(8) 



The only constraint which fixes the coefficients in front of the WZ counterterms is gauge invariance. 
Specifically, the anomalous variation of San is compensated by the variation of Swz- Imposing this 
condition one discovers that the scale of the WZ counterterms (M) becomes the Stiickelberg mass term 
Mst = Ml. This is found in the defining phase of the model, in which the realization of the gauge 
symmetry is in the Stiickelberg form. Obviously, in this phase only the B gauge boson is massive (in a 
Stiickelberg phase). The breaking of the electroweak symmetry, triggered by the Higgs potential and 
the transition to the mass eigenstates determines a rotation of the Stiickelberg axion b into a physical 
axion x plus some Nambu-Goldstone modes. This rotation brings in a redefinition of the suppression 
scale M, which now coincides with the mass of the extra Z' gauge boson, as shown in an appendix. 

2.2 The scalar potentials and their axion-dependent phases 

In previous studies it has been shown that anomalous abelian models, realized in the case of poten- 
tials with 2 Higgs doublets, both in the non-supersymmetric and in the supersymmetric cases, are 
characterized by the presence of an axion-like particle in the spectrum. In the context of the 2 Higgs 
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doublets model shown in detail in [36\ [38] the presence of PQ-breaking terms in the scalar potential 
allows the axion to become massive. The PQ symmetric contribution is given by 

VpQiH^, Hd) = ^ (MHa + laaiHiHaf) - 2lud{HiHu){H\Ha) + lUHlr^HdW (9) 

a=u,d 

which is a pure Higgs scalar potential, while in the PQ-breaking terms we introduce a dependence on 
the axion field h by means of explicit phases 



Vf. (5 {Hu,Hd,h) = bi [HlHae «i + [HlHde 



+A2 (hIHu) (HiHae ) + A3 (H^Ha) (HlHae '^h) + h.c. 



(10) 

where Aq^ = 1u ~ Qd ^ ^1 has mass squared dimension, while A^, A2, A3 are dimensionless couplings. 
In the scalar potential we can isolate three sectors, namely, two neutral and one charged sector, which 
are described by the quadratic expansion of the potential around its minimum 



VcP-eveniHu, Hd) + VcP-odd{Hu, Hd, b) + V±{Hu, Hd) 



(11) 




+ {lmHj'jmHd\ai)M^ ImHd'^ \. (12) 



• The Charged Sector 

In the charged sector we find a zero eigenvalue of the mass matrix, corresponding to the Goldstone 
mode and the nonzero eigenvalue 

mjj+ = 4X'udV^-2( '^.^ +2h+tanpi2 + cotPh)v'^, (13) 

\u^sm2p / 

corresponding to the charged Higgs mass. The two vevs of the Higgs sector are defined by Vd = 
V cos P]Vu = v sin /3, with v"^ = v"^ + v"^. The rotation matrix into the physical eigenstates is 



Hu+ \ _ fsinp -cos/?\ / 
Hd"" ) ~ Vcos/3 sin/? J \ 



(14) 



The CP-even Sector 
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In the neutral sector both a CP-even and a CP-odd subsectors are present. The CP-even sector is 
described by N2 which can be diagonahzed by an appropriate rotation matrix in terms of CP-even 
mass eigenstates {h^,H^) as 



with 



and 



tan a 



sm a — cos a 
cos Q sin a 



AA2(1,1)-AA2(2,2)-^ 
2AA2(1,2) 



(15) 

(16) 
(17) 



A = (AA2(1, 1))2 - 2AA2(2, 2)AA2(1, 1) + 4 (AA2(1, 2))^ + (AA2(2, 2))^ 

The definition of these matrix elements is left to an appendix. The eigenvalues corresponding to the 
physical neutral Higgs fields are given by 

1 



mlo 



AA2(1,1)+AA2(2,2)-VA 



AA2(1,1)+AA2(2,2)+VA 



(18) 



We refer to [36] for a more detailed discussion of the scalar sector of the model with more than one 
extra U{\). 

• The CP-odd sector 

The symmetric matrix describing the mixing of the CP-odd Higgs sector with the axion field b is 
given by A/s. After the diagonalization we can construct the orthogonal matrix that rotates 
the Stiickelberg field and the CP-odd phases of the two Higgs doublets into the mass eigenstates 
(X,G0,G0) 



(19) 



The mass matrix of this sector exhibits two zero eigenvalues corresponding to the Goldstone modes 
G]*, G2 and a mass eigenvalue, that corresponds to the physical axion field Xi with a value 



/ /mi/0 ^ 




( ^ \ 




= 0^ 


G? 









1 + 



^ V sin 2/? 



with the coefficient 



A/i 



4 4Ai + A3 cot P + 



1 + 



Ml 



hi 2 
v"^ sin 2(3 



+ A2 tan/3 



(20) 



(21) 
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The mass of this state is positive if < 0. The Goldstone bosons {Gz, Gz') are obtained by orthonor- 
mahzing that span a two dimensional space. Notice that, in general, the mass of the axi-Higgs 

is the result of two effects: the presence of the Higgs vevs and the presence of the Stiickelberg mass via 
the PQ-breaking potential. In the particular case of a charge assignment such that = q^, in the PQ- 
breaking potential the dependence on the axion field disappears {Vf> ^ Hd, b) ^ Vf> {Hu, Hd)) 
and the rotation matrix simplifies to 



(22) 



For this particular assignment of the Higgs charges the Z and Z' bosons are still massive, as can be 
seen from eqs. (jl82| I183p . A brief counting of the physical degrees of freedom shows, also in this 
case, that we expect only one physical particle in the CP-odd sector. Then, in this particular case, 
it is easily found that the model doesn't exhibit Higgs-axion mixing because the physical degree of 
freedom , as identified by the scalar potential, is a combination of the imaginary parts of the two 
Higgs /mi7°, ImH^, while the axion is only part of the Goldstones modes Gz and Gz', identified by 
an inspection of the derivative couplings. 




— COS (3 


sin/3 




/ \ 


sin (3 


cos j3 


:) 


G\ 













3 Axions from the decoupling of a chiral fermion 

Other realizations of these effective models are obtained by studying the decoupling of a chiral fermion 
from an original anomaly- free theory, due to large Yukawa couplings [ID]. The remnant axion, in this 
particular realization, is the surviving massless phase of a heavy Higgs. We will illustrate briefly this 
approach sketching the derivation, though in the case of a simple model, in a section below. Obviously, 
in these types of completions of the anomalous theory, the challenge of the construction would consist 
in the identification of a pattern of sequential breaking of the underlying anomaly-free theory in order 
to generate suitable axion-like Wess-Zumino interactions, which are not part of our simple example. 

For instance, considerable motivations for this reasoning comes from unified models based on an 
anomaly-free fermion spectrum assigned to special representations of the gauge symmetry. Specifically, 
one could consider the 16 of 50(10) in which find accommodation the fermions of an entire generation 
of the Standard Model plus a right handed neutrino. The decoupling of a right handed neutrino could 
leave a remnant pseudoscalar in the spectrum with axion-like couplings. While the explicit realization 
of this construction and the (sequential) breaking of the original GUT towards the spectrum of the 
Standard Model is rather complex, the implications of these assumptions can be grasped by a simple 
model. 

To illustrate these points, we introduce a simple toy model and show step by step that a specific 
form of the decoupling can generate a certain dynamics at low energy which is completely described 
by an effective action with Stiickelberg and a Higgs-Stuckelberg phases, Wess Zumino interactions 
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and higher dimensional operators suppressed by the Stiickelberg mass. It should be mentioned that 
in our example, the low energy gauge boson B, which has anomalous effective interactions, would be 
massive in the Stiickelberg form. We recall that the study of the Stiickelberg construction has been 
discussed recently in several works |4HI42j (see also |43j) for non-anomalous theories, with its possible 
experimental signatures. 

The model requires two Higgs fields, here assumed to be two complex scalars, and a potential 
characterized by a first breaking of the anomaly- free gauge symmetry at a certain scale (v^f,), followed 
by a second breaking at a lower scale vh {vh « v^). The heavy Higgs is assumed to decouple 
(partially) after the first breaking. Specifically, the decoupling involves the radial fiuctuations (p) of 
the field (/>, and all the interactions which are characterized by operators which are suppressed by a 
certain power of p/v^. We expand the heavy Higgs as 



with 9 denoting a massless phase that may be rendered massive during the process of decoupling of 
the radial excitation by some small tilting, as it occurs for the ordinary Peccei-Quinn axion (PQ). The 
(almost massless) phase remains in the low energy theory. The Stiickelberg axion is identified from 
^ in a certain way, that will be specified below. Also we assume, for simplicity, that only one chiral 
fermion becomes heavy in the course of decoupling of the heavy Higgs, and is integrated out of the 
low energy spectrum. As we have already stressed, our approach can be made more realistic, but we 
expect that the crucial steps that bring to its specific effective action at low energy can be part of a 
more complete theory. 

The Yukawa couplings, expanded around the vacuum of the heavy Higgs, show the presence of a 
complex phase {6) that we try to remove by a chiral redefinition of the integration measure before we 
integrate out the heavy fermion. It is this chiral redefinition of the fermionic measure which induces, 
by Fujikawa's approach, typical Wess-Zumino terms in the low energy effective theory. This theory, 
obviously, admits a derivative expansion in terms of the large scale v^, which can be systematically 
captured by a derivative expansion in l/f<^, or equivalently, the Stiickelberg mass, since the two scales 
are related (Mi gs^tf,). 

3.1 Partial integration 

To be specific, we consider a model with 2 fermions and a gauge symmetry of the form U {1)a xU{1)b, 
where A is vector like and B is the anomalous gauge boson. We define the lagrangean 




(23) 




i=l 




(24) 
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Uil)B 
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Iar 
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(2) 
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(2) 

Iar 


(2) 

Tbr 
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1a 


Qb 





<t> 
lA 


fB 



Table 1: Charge assignments for the A-B toy modeL 



('2) (2'] 

where we have neglected the Yukawa couphng of the hght fermion(s) V'}, , > which are proportional 
to the vev of the light Higgs vh- For simplicity we may consider a simple scalar potential function of 
the two Higgs (f) and such as V{4),H), that as we have mentioned, admits vacua which are widely 
separated. While this would induce a hierarchy between the two vevs, and could be the real difficulty 
in the realization of this scenario, one possible way out would be to consider V{(j), H) to be the sum 
of two separate potentials. Since the phase of the heavy Higgs survives in the low energy theory as a 
pseudo-goldstone mode, it may acquire a mass if the potential in which it appears is tilted. 
We show in Tab. [T]the charge assignments of the model. We define 

Df.i^? = {^^, + iq%gAA^ + iqBL9BB^^^p]/;\ (25) 

(26) 

Under a gauge transformation we have ^/^ — > V' 

^J) ^ e-"^n9Be^(^) (27) 

with 6B^ = B'^-B^ = -d^e. 

We assume that the charge assignments are such that the model is anomaly-free. Notice also that 
B, in this realization, becomes massive via a first breaking at the large scale and then its mass gets 
corrected by the second breaking, characterized by the scale vh- 

We parameterize the fluctuations of the field (p around the first vacuum in the form 



^Jt±Pe-<<^Be (28) 
v2 
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from which we obtain the first contribution to the mass of the B gauge boson in the form Mi = q'^gBV^- 
As we are going to show next, this mass can be taken to be the Stiickelberg mass of a reduced Higgs 
system if we neglect the radial excitations. In fact we have 

\D,<t>\' = \ {d.cpf + (^^) ' {qUb? {-d,e + B,f , (29) 

and we isolate from the phase 9 of this exact relation a dimensionful field b which will be taking the 
role of a Stiickelberg mass term as 

& = (30) 



qB9BVcf> 



We can expand ()29p in the form 



1 



\D^<P\^ = ^id>^- M^B^f + 0{p/v), (31) 

with Ml = q'^gBV<j>, defined to be the Stiickelberg mass. The decoupling of the radial excitations of 
the very heavy Higgs from the low energy lagrangean generates a Stiickelberg mass term on the rhs 
of (j3ip . whose phase 9 is at this stage massless. Notice that after the second symmetry breaking, the 
mass of the B gauge boson will acquire an additional contribution proportional to QeQe'^h, in analogy 
to the first breaking, that is 



Mb = ^Ml + [gsq'^VHf. (32) 

Notice also that after the first radial decoupling of the heavy Higgs the Yukawa mass terms are 
affected by a phase dependence that can be eliminated from the effective lagrangean via an anomalous 
transformation. To illustrate this point consider the expansion of the Yukawa term around the vacuum 
of the heavy Higgs 

AV^i'VV'g^ = ^^{v^ + p)i^^l^i^^i^e-<3Be (33) 

which is affected by a phase that we will try to remove in the course of the elimination of the heavy 
degrees of freedom of the mother theory. Notice that in this case we do not take a large Yukawa 
coupling (A), as in previous analysis |441 [15] . since the large fermion mass of ■0*-^^ is instead obtained 
via the large vev of the heavy Higgs, v^. For this reason, having defined the Stiickelberg mass Mi in 
terms of the same vev, after neglecting the radial contributions we obtain 

A^^VV'SJ^ = KMiV5i'V2^e-''B9Be, ^ = ^q^QB. (34) 

Before performing the partial integration on the heavy fermion , it is convenient to define a change 
of variables in the functional integral, in order to remove the phase-dependence on 9 present in the 
Yukawa couplings. For this reason, let's consider the part of the partition function directly related to 
the heavy fermion which is involved in the procedure of partial integration. This is given by 
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where 



Z^^\A, B) = j V^iWiWiWu^e^^ (35) 
£(1) = V'W^ V'?) + Vig^^ + KMiV^i^VSJ^e-^^sffl^^ + h.c. (36) 



and we have neglected the contributions proportional to the radial excitation of the heavy Higgs. At 
this point we try to remove the phase 9 from the Yukawa couplings by performing a field redefinition 
in the functional integral of the heavy fermion. We set 



41 = e--^^i-VS 

^g) ^ e--^^l.-VS, (37) 

where from gauge invariance we have 

C + 4 - Q^bI = 0- (38) 
The field redefinition induces in the integration measures two jacobeans 

V4^V4^ = JrD4^v4r~^ (39) 
which are computed using Fujikawa's approach (see for instance |46] ) . We obtain 



= e-^'?^H^<^^^^)«. (40) 



In this case F^^l^r = [D^, D^]^ contains both gauge fields {A, B) and the corresponding gauge 
charges of the heavy (L, R) fermions such as, for instance, 

F,uL,R = iq%RF^u + iqilRF^i'u- (41) 
The structure of the effective action after the field redefinition takes the form 

Z^^\A,B) = J v4^v4^v4^V,l;f^e'I'''^''"''+^^- (42) 

where 

C"^'^ = 4^ - ^^21^ O) 4^ + i,f [lj> - e) 4^ + nM.^^i^f + h.c. (43) 

with the Wess-Zumino (WZ) lagrangean obtained from the expansion of the OF A F terms. These 
are suppressed by the Stiickelberg mass term Mi {6 = h/Mi). 
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At this point we can perform the Grassmann integration over the heavy fermion, which trivially 
gives the functional determinant of an operator, P, exphcitly given by 



/ 



(44) 



where v'^ = v'^ /^/2. The remaining terms in the total partition function of the model can be obtained 
from the functional integral 



where 



with 



C'^^^ + Cwz + Tr log P + ^ (a^ - MiB^f + \D,,H\'' - V{H, 6) (46) 



^'^'^ = -\fI - \fI + ^ + ^g^^ (47) 

The derivative expansion of the effective action can be organized in terms of corrections in the 
Stiickelberg mass. Obviously, a similar approach can be followed for the integration of a Majorana 
fermion, which is slightly more involved. The basic physical principle, however, remains the same also 
in this second variant. In this case the functional determinant can be organized as in |47j . 

There are some implications concerning the two realizations of this class of effective actions, es- 
pecially in regard to the possible mass of the axion as a dark matter candidate in the various models 
that share the effective actions that we have presented. The first observation concerns the absence of 
a direct Yukawa coupling between the heavy Higgs and the light fermion spectrum, which is part of 
the effective action after partial integration on the heavy fermion modes. This feature is absent in the 
MLSOM, and turns out to be rather important since it affects drastically the lifetime of the axion, as 
we are going to elaborate in the following sections. We will find that a GeV axion is favoured by the 
mechanism of partial decoupling but is not allowed in the MLSOM. In this second case a very light 
axion is necessary in order to have a state which is long lived and that can be a good dark matter 
candidate. 



3.2 Parametric solutions of the anomaly equations 

It is clear that the typical effective action isolated by the decoupling of (one or more) chiral fermions 
can be organized in terms of the defining lagrangean plus the WZ counterterms, which restore the 
gauge invariance of the model. Therefore, up to operators of mass dimension 5, the two lagrangeans are 
quite overlapping at operatorial level. For this reason, we will construct a complete charge assignments 
for these models, starting from an anomaly-free theory, with a spectrum that we deliberately choose 
to include one right-handed neutrino per generation, and which we will decouple from the low energy 
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dynamics according to the procedure described above. Of course, other choices are also possible. As 
we have already stressed, the motivations for selecting this approach are not just of practical nature, 
although it allows to generate effective anomalous models with ease. For instance, one could envision 
a scenario, inspired by leptogenesys, which could offer a realization of this decoupling mechanism, 
although its details remain, at the moment, rather general. We will not pursue the analysis of this 
point any further, and leave it as an interesting possibility for future studies. However, we will discover, 
by using the decoupling approach, that a significant class of charge assignments of intersecting brane 
models can be easily reproduced by the free gauge charges which parametrize the violation of the 
conditions of cancellation of the anomaly equations. We should also mention that the dependence of 
our results on the various charge assignments is truly small, showing that the relevant parameters of 
the models are the Stiickelberg mass, the anomalous coupling and the parameters of the potential, 
which control the axion mass in each realization. 

To proceed, we impose first the conditions of cancellation of the gauge and of the mixed gravitational- 
U{1)b anomalies, thereby fixing the U{1)b charges, followed by the conditions of invariance of the 
Yukawa couplings, in order to determine the charges of the two Higgs |48] . We take the 17(1)^ fermion 
charges to be family-independent in order to avoid possible constraints from flavor-changing neutral 
current processes. We label the generic fermion charges under the additional group U{1)b as shown 
in Table [3 



Ql 


UR 


dR 


L 




^R 








ql 




q^ 

qu 



Table 2: Labels for the gauge charges of the fermion spectrum. 

For every anomalous triangle we allow, in general, a WZ counterterm whose coefficient has to be 
tuned in order to satisfy the conditions for anomaly cancellation. For the fermion charges , 
we find the following constraints 

BSU{2)SU{2) : + 3g|^ - Cbww = 0, 

BSUi3)SUi3) : + < - 2q^^ - CBgg = 0, 

BYY : 3gf^ + Qq^^ + Sg^ - ^Cbww - Cbyy - CBgg, (48) 

where the coefficients appearing in front of the WZ counterterms are proportional to the charge 
asymmetries 

Cbww oc ^ 
/ 

^Bgg (X^e^, 
Q 
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Cbyy oc ^ ^/^^) 

/ 

CBBB^Y^Of^^, (49) 

/ 

which are detailed in an appendix, and with the hypercharges of U{1)y given in Tab. (I3D- 

If we consider the charges qQ^,qf as free parameters of the model, Cbww i Csgg , Cbyy can be in 
principle expressed in terms of these parameters. The other three conditions coming from the gauge 
invariance give the following further constraints 



B \2 of„B \2 , Q^„B^2 Q/„B \2 , c^„B ^2 

BBB 

9ql + 9g^^ + 3gf^ - 6gf - 18q§^ - 3Cbgg = 0, (50) 



YBB : -3{qir - 3(gfJ^ + 3{qtr - 3(g0J^ + 6{qSJ' - Cybb = 
%Qlf + 3(gfj' - 6(gf )' - 18(g|j3 + 9(gfj3 _ ^^^^ ^ q 



where the condition on the BRR triangle comes from the mixed gravitational-C/(l)B anomaly cancel- 
lation. From the gauge invariance of the Yukawa couplings (see Lagrangian (j73p ). we obtain 





9d 
2 


< 


= 0, 




Qu 
2 


< 


= 0, 




Id 
2 




= 0, 




d + 


q^ 
2 


= 0, 



(51) 

which can be used to constrain the charges of the two Higgs doublets qu^q^ ^"^^ ttis counterterms 
Cbww iCbqq- Collecting the constraints in eqs. ([5T]) . (|18]) and ([50]) we obtain a set of ten equations 
whose solution allows us to identify a class of charge assignments that we call / 

/(g|^,gf , Ag^) = (g|„<;<,Qf ,gf«,gf ,9^ )• (52) 

These depend only upon the three free parameters qq^, q^, Ag^, where Ag^ = qS ~ Qd ■ The explicit 
dependences are shown in Table [3l while the related WZ counterterms take the form 

CBYY = -l{qE- 5q^, ) + 2 Ag^ , (53) 

Cybb = Kd? - ^ [l8(g|j^ + 8gg,A,^ + (A,^)^] , (54) 

Cbbb = -6{qEf + n{q^,f + 72{q^j'Aq^ + 18g|jAg^)2 + ^iAq^f, (55) 

CBgg = \Aq^, (56) 
CBWW = qE + 3q^^, (57) 
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f 

J 






U(1)y 




Ql 


3 


2 


1/6 


R 


ur 


3 


1 


2/3 





dn 


3 


1 


-1/3 




L 


1 


2 


-1/2 




eR 


1 


1 


-1 




^R 


1 


1 








Hu 


1 


2 


1/2 






1 


2 


1/2 





Table 3: The three-parameter family /{qq^, qj^ , Ag^) of solutions for fermion and scalar charges. 

where in particular, from the charge assignment shown in Tabled we identify the counterterm for the 
mixed gravitational-C/(l)B anomaly with 

CBGG = 2{-qE+q^^+Aq''). (58) 

Then the WZ counterterms, as defined in general in eqs. (jl94p . can now be specialized in terms 
of the different charge assignments f{qQ^,q^^,Aq^), just by substituting the corresponding chiral 
asymmetries. This function will appear in several of our plots. 

Finally, since in the case q^ —q^ =0 the matrix would become trivial, we require the following 
relation between the Higgs charges 

Qu-Qdy^O (59) 

where, in particular, q^ — qj^ = A is exactly the value implied by the charge assignment derived from 
the Madrid Model (see Table [6|) for the two Higgs. We will be using this value to constrain the chiral 
asymmetry 0^ by means of eq. (j58p . and will be taken as the starting value for all our comparisons. 
Notice that the family f{qQ^,q^,Aq^) for the particular choice qq^ = —l,qi = —1 reproduces the 
entire charge assignment of the Madrid Model 

/(-1,-1,4) = (-1,0,0,-1,0, +2, -2). (60) 
3.3 The Madrid model 

We just recall, as already mentioned, that the charge assignment for our anomalous (brane) model 
that we consider is obtained from the intersection of 4 branes (a, b, c, d) with generators {qa, qt, qc, qd) 
which are rotated on the hypercharge basis U{l)xi with i = A,B,C and U{1)y, with an anomaly 
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free hypercharge. The U{l)a and U{l)d symmetries are proportional to the baryon number and the 
lepton number respectively. The U (l)c symmetry can be identified as the third component of the 
right-handed weak isospin, while the U{l)b is a PQ-like symmetry. A detailed discussion of this 
construction can be found in [4J and [49j. The identification of the generators involve the solution 
of some constraint equations. In general, for a simple compactification the solutions of these 
equations are parametrized by a phase e = ±1, the Neveu-Schwarz background on the first two tori 
f3i = 1 — bi = 1, 1/2, the four integers na2,nhi,nci,nd2 which are the wrapping numbers of the branes 
around the extra (toroidal) manifolds of the compactification, and finally a parameter p = 1,1/3. 
One of the possible choices for these parameters is reported in Table H] which identifies a particular 
class of models, the so called Class A models. The result of this D-brane construction is the charge 



V 






na2 




rici 


nd2 


1/3 


1/2 


1 


na2 


-1 


1 


1 - na2 



Table 4: Parameters for a Class A model with a D6-brane . 



/ 


Ql 


UR 


dR 


L 






Qy 


1/6 


2/3 


-1/3 


-1/2 


- 1 





1b 


-1 








-1 









Table 5: Fcrmion spectrum charges in the F-basis for the Madrid model [49] . 

assignment specified in Table [5] whose corresponding fermion spectrum is anomalous under the extra 
U{1)b abelian symmetry. Imposing the gauge invariance of the Yukawa couplings, see eq. (j73p . we 
constraint the charges of the Higgs doublets to the values specified in Table [H 





Y 








1/2 





2 


Hd 


1/2 





-2 



Table 6: Higgs charges in the Madrid model. 
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4 Trilinear and quadrilinear interactions of the axi-Higgs from the 
MLSOM scalar potential 



One of the objectives of this work is to quantify the decay rates in the various channels of the axi-Higgs 
X and of the two Higgs bosons and of the CP-even sector, and to explore some possible channels 
in which the production of an axi-Higgs can be realized at the LHC. For this goal we proceed with 
a careful inspection of the interaction lagrangian, in order to extrapolate all the relevant couplings 
and interactions of the axi-Higgs and of the CP-even sector with the other particles. We start this 
analysis by collecting first all the trilinear and quadrilinear interactions of the axi-Higgs that emerge 
from the scalar potential and then move to the mixed vertices which involve both the CP-even and 
CP-odd sectors. 

Collecting the quadrilinear vertices we obtain 



X 



(61) 



where we have defined 



R^ 



i?4 



'2Mi 



X \31 



(62) 



The first contribution (i?i) is extracted from the diagonal part of the Higgs potential (i.e ~ laaiH^Haf), 
the second originates from the non-diagonal u-d terms (~ lud{HuHu){H^^Hd)), the third comes from 
the contribution of the PQ-breaking potential proportional to li , while R^ is the contribution of the 
last two pieces of the same potential which are proportional to I2 and ^3. 

The quadrilinear couplings of the axi-Higgs with the neutral Higgs sector involve interactions 
between two axions and the two neutral states {H^,hP). We can write the interaction lagrangian as 
follows 



+ 



R 



X^HOHO j^x^HORO j^x^H'^HO j^x^H^H^ 
+ K2 + -ft3 + -K^ 



R 



+ Ri 



+ R1 



(63) 



where the coefficients R^ are defined in an appendix. 
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The trilinear interactions of the axi-Higgs with the neutral Higgs sector exhibit couphngs with two 
axions and one Higgs state H'^, hP. The interaction lagrangian can be written as 



"-X^higgs = ^x^H" + ^x'^h" 



(64) 



where we have defined 



jC.y2fjO 



.1=1 



i=l 



2/1° 



2 1.0 



(65) 



Again, the Rf ' coefficients are hsted in an appendix. It is important to note that these couphngs 
are also present in a general 2HDM, while they are absent in the MSSM due to the strong constraints 
obtained by imposing supersymmetry. 

4.1 Self interactions in the CP-even sector 

The self interactions of Hq and fiQ can be described as above, by analyzing the quadrilinear and 
trilinear vertices generated by the rotation of the fields in the physical basis after electroweak symmetry 
breaking (EWSB). Starting from the quadrilinear interactions we can write 



(66) 



where 



tt4 tt4 tt4 zt4 






h2zj2 !,2rj-2 !,2u-2 

<° ° + <° ° + Ry ° 



R 



R 



Hoh^ 



+ R, 



■3 

0^0 



(67) 



tt4 

The coefficients Rf can be found in an appendix. Also here it is interesting to observe that Ri and 
i?2 are in general related to the PQ symmetric part of the scalar potential, while Rz and i?4 come 
from the PQ-breaking terms. 

The trilinear interaction lagrangian can be written as 



ho' 



-ho^Ho 



+ c 



(68) 



where we have defined 



tt3 tjZ tjZ tj3 ' 



Hi 
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i3 ^3 i3 i3 

rI" + R2" + R3" + R2" 



R 



hlHo 



+ R. 



hlHo 



R 



+ R. 



(69) 



All the coefficients Rf ^ are given in an appendix. 



4.2 Trilinear interactions of the CP-even sector with the and Z gauge bosons 

Since, in general, the branching ratios for the decay of the Higgs into a pair of vector bosons 
or ZZ are relevant in a certain kincmatical regime, it is important to quantify the tree level decay 
rate for this channel, and to give an estimate of the coefficients of the trilinear interactions of Hq and 
ho with two gauge bosons W^W^ and ZZ. For the charged it is straighforward to obtain the 
corresponding coefficients 

Cww = Y (sin (X Vd- cos a Vy) , 
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Cww = Y (sin a ■"d + cos a Vy) 



(70) 



The calculation of the coefficients for the analogous interactions with the Z's is more complicated 
because of the structure of the model. For this purpose it is useful to introduce the following coefficients 



/i = 



fl + h (^/i - 2gBqixB) + 2xb [xsg^ + QBqi (gBq%XB - ^J fl + 4.g^xl) 



6 = 



/i + /i (y7F+V4 - '^dBQlxB) + 2xB [xBQ^ + QBql [gBqlxB - ,/f[+4^) 



2V2(Ag^xl + h^fl + Ag^xl) 



and the interactions H-Z-Z at tree level - summarized by the coefficients C^z - a^re given by 

Czz = {vd sin a-Vug'^ ^1 cos a) , 
Czz = g^ il sin a + Vdg^ cos a) , 



(71) 



(72) 



where 5^ = + 52- 
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5 The Yukawa couplings and the axi-Higgs 



The couplings of the two Higgs and of the axi-Higgs to the fermion sector are entirely described by 
the Yukawa lagrangian. The Yukawa couplings of the model are given by 

C^rut. ^ -T''Qj^HddR-V^dRH\QL-T^Qdi^2K)uR-T^UR{ia2H*jQL 
-T'LHaeR - T^eRH\L - V'L{za2H*JuR - T^UR{ia2H*jL 
= -V^ dH^PRd - r'^ dH^/Pid - T" uH^*Pru - uH^Plu 

-r^ eHl>PRe - eH^/PLe - VHI*Pru - T'' VHI^Pli^, (73) 

where the Yukawa coupling constants V^, F", F*^ and run over the three generations, i.e. u = {u, c, t}, 
d = {d,s,b}, u = {ue, f^, i^t} and e = {e, /j,, r}. Rotating the CP-odd and CP-even neutral sectors 
into the mass eigenstates and expanding around the vacuum we obtain 

Hu = vu + ^ 

^ ^ (/iOsina-gOcosa)+z {Of,X + O^G? + O^sG^) ^^^^ 

^ {h^ cos a + H^ sin a) + i {0^,x + 02^2^1° + 0^3^°) ^^^^ 

so that in the unitary gauge we obtain 

= Vu + ^[{h" sin a- cos a) +iO^-^x] 

= Vu + ^ [{h° sin a- cos a) -iN cos (i x] (76) 



V2 

= Vd + ^ [(/i° cos a + sin a) + i O^^ x] 
V2 



Vu + ^ [{h° cos a + sin a) + iNsmP x] , (77) 



V2 

where the vevs of the two neutral Higgs bosons Vu = v sin /3 and = v cos (5 satisfy 



tan/3 = ^, v = Jvl + vl (78) 

Vd ^ 

We have also relied on the definitions of introduced in a previous work [38] 

0\^ = -N cos 13, (79) 



23 



0^^=Nsmf3, (80) 

that we have reported in an appendix. For convenience we have introduced the fohowing normahzation 
coefficient 



A^= , (81) 



The fermion masses are given by 



rU T^U 

md = VdT'^, rUe = VdT'', (82) 

where the generation index has been suppressed for brevity. The fermion masses, defined in terms of the 
two expectation values Vu, Vd of the model, show an enhancement of the down- type Yukawa coupHngs 
for large values of tan j3 while at the same time the up-type Yukawa couplings get a suppression. The 
couplings of the hP boson to fermions are given by 



/ sin a , n\ , . 

-r''uLUni^-j=-h'^j +C.C. (83) 

The couplings of the boson to the fermions are 

0\ _ T-d-r r I sina q\ „. I cosa o\ T^e^„ I ^^^oi q 



-r^I7..^(-^if°)+c.c. (84) 
For later reference we group together the couplings of the axi-Higgs x with the fermion sector 

^Ynk[X) = -r dLdR ^ X j - r ULUR ^ X j - r ClCR ^ X 

^u- ( . NCOS (3 \ 
-T VLVR y-i—j^x j + c.c. (85) 

We have listed these couplings in Tab. d?]) where the normalization coefficient N is defined in (jSip . 
From the Yukawa couplings of eq. (173]) and relations (f71|l . ([751) we can extract the coupling of the 
Goldstone boson to the fermions 

/:Yuk(G2°) = -T^dl^^Gl^PRd-T'^-dl^-i^Gi^PLd-T-ul^-i^Gi^^ 

-r-n p,u - e PRC - P,e 

-r^I7 Pr. - T^v(^^Gi) P,u. (86) 
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up-fermion 


down-fermion 


Higgs SM 


V 




Lighter Higgs /i° 


— sina/smp 


— —^ COS a/ cos p 


Heavier Higgs 


^ cos a/ sin/3 


— ^ sina/ cos/3 


axi-Higgs x 




-z^ ATtan/? 



Table 7: Couplings of the neutral MLSOM Higgs bosons to up- and down-type fermions, and comparison with 
the fermion couplings of the SM Higgs boson. 



Using the expression of we can compute the coupling between the Goldstone boson and the 
down-like quarks that takes the form 



- r'* dR^PR d-V^ dH^*PL d 



V2 



-N 



Ml 



id-f^'dG^ 



(87) 



and similarly for the other generations. These expressions have been used in order to fix the explicit 
form of the Wess-Zumino (WZ) counterterms using the condition of gauge invariance. 



6 Decay rates of the axi-Higgs 

We proceed to compute the partial decay widths and the branching ratios of the axi-Higgs for different 
decay modes in the CP-odd sector of the MLSOM, taking the mass of the axion as a free parameter. 
As we have already mentioned, in the case of the MLSOM, there is an interesting window in which 
the axion acquires a lifetime typical of a good dark matter candidate. This mass value, which is the 
same as that of a traditional Peccei-Quinn axion 10^^ eV, or in the ultralight mass window), is 
not the most interesting one for studies of this particle at the LHC. The reason of this result has to be 
found in the fact that the most relevant channels for the production of a particle of this mass are 1) 
the pscudoscalar vertex with a top or bottom quark loop (the dominance of one or the other fermion 
contribution depends closely on the value of tan/3); 2) the direct WZ vertex in which the axion is 
radiated off by a gauge field. The WZ term is quite small compared to the contribution from the 
fermion loop, which is instead dependent on the mass of the axion. For an ultralight axion the loop 
contribution is rather small and the chances of producing a particle of such a mass by gluon fusion or 
in qq annihilation of light quarks are quite small. For this reason, if we are interested in the study of 
a GeV axion, which is the goal of the numerical sections that follow, we are automatically excluding 
a long-lived particle. On the other end, in this mass region, we are instead analyzing a particle whose 
behaviour is Higgs-like but with a direct (although small) direct coupling to the gauge fields. At the 
same time, the Higgs-like nature of the axion can be investigated by taking its mass in the several 
GeV region, say in the 100-120 GcV range. Our results, however, are quite general, in this respect, 
and can be used for direct studies of this particle in any mass range. As we have already stressed. 
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what makes a distinction between a "standard" CP-odd Higgs state and the axion of the MLSOM are 
the WZ interactions, which are, in any case, subdominant compared to the triangle diagram in any 
mass range. 

In the case of fermion decouphng one can proceed with similar considerations, although the con- 
clusions are rather different and will be addressed below. We will describe in a final section the main 
properties of the axion if its origin is to be traced back to a decoupled Higgs sector, which show, in 
this second realization, that the axion can be long lived and with a mass in the GeV range. 

The relevant parameters which appear in the decay are the following coefficients 



Vu V tan p 



-I — - 


qX _ 


-i—!^NtanP, 


^d 


V 


I— 


QX _ 


in 

-i— Ntanf], 




V 



V tan (5 

which will be essential in order to establish the size of the various decay channels. 

Since we are interested in a relatively light axi-Higgs, we have focused our study on a kinematical 
mass range going from 1 to 100 GeV. The fermionic decay channels that we consider are the 66, cc, 
ss for the tree level decays into quarks, rf and ^Jl for the decays into leptons. At one-loop order we 
consider the decay into two photons, two gluons and in one photon and one Z boson. We have added 
both the massless contribution coming from the WZ counterterm and the fermion loop contribution 
from a pseudoscalar triangle. The total decay rate of the axi-higgs in this approximation is given by 



X 

ir 

q=s,c,b l=IJ',T 



• The tree level decays into fermions: x ~^ ff 

At leading order, for the tree-level process x ~^ ffi we obtain the decay rate 



n^^ffl^'^e'Q}ic^-'f-N,U)^l-[^)\ (90) 

for a value of the fermion mass below the pair production threshold (Atjij < m'^). The pseudoscalar 
couplings to the fermions (c-^'-^) have been defined in Eq. (j88p . 

The leading decay is x ~^ bb, due to the suppression of the fermion couplings of the up-type 
fermions (clearly shown in Table [6]). We show the variations of the branching ratio (BR) of the 
pseudoscalar for different charge assignments /(— 1, — 1, Ag'^), and as observed before, there are no 
substantial differences induced by the selection of different assignments. 

• One-loop decays into photons and gluons: x ^ 77 X ^ 99 
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We now compute the partial decay width of the axi-Higgs boson into two photons x ~^ 77- The 
invariant matrix element considered for the process is the sum of the two contributions shown in 
Fig. [TJ The first amplitude (Fig. [1^) is a massless WZ vertex 



(91) 



where the coefficient g^^ comes from the counterterm given in formula (jl92p . The second amplitude 
(Fig. [Da) is a pure massive contribution 

M'}''{x^ll)=J2N,{f)iCo{ml,mf)c^:ls[fi,u,h,k2], f = {u,d,u,e} (92) 
/ 

where Ndf) is the color factor, 1 for leptons and 3 for quarks. In the domain < < 2m f the 
pseudoscalar triangle when both photons are on mass- shell kl = k^ = is given by the expression 



Co(m^,m/) = 



mt 2 
arctan 



nif 2 
arctan 



2m i 



(93) 



with 



2m 



m^ 



while in the domain 2m f < m^. it becomes 

Co{m^,mf) = ReCo(m^, m/) + iImCo(m^, mj). 



Here we have set 



ReCo(m^,mj) 



ImCo(m^,m/) 
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In the numerical analysis presented below, we have introduced the function /(r), defined in any 
kinematic domain, whose real part is given by 



Re[/(T)] = I 
while its imaginary part is 

Im[/(r)] 



(arcsin l/-yr)^ 

^2 f l + v^l-T 



vr 



ifr > 1 
ifr < 1 



(98) 



ifr > 1 
ifr < 1 



(99) 
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where r = Arnj / . 

X Finally, the 1-loop decay x — > 77 is given by the following amplitudes 

A^^'^(x-77)=A^lV'z + -^r 
and the rate computed from the two contributions shown in Fig. [1] is 



(100) 



mr 



r(x-77) = ^< 



8(9?-,)' + 5 



(101) 



In Fig. [T^ we have isolated the massless contribution to the decay rate coming from the WZ counterterm 
xF^F^ whose expression is 

,3 

(102) 




(o) (b) 



Figure 1: Massless plus massive contributions to the x ^ 11 process. 

We should notice that the massive contribution from amplitude ()92p is completely independent 
of the anomalous coupling gs, which does not appear in the coefficients c-^'-^, as can be seen from 
Eq. (j88p . For the decay into two gluons we proceed in a similar manner (see Fig. [8]) and the amplitude 
is given by 

M'^ziaa ^ x) = "^g^gsifi, u, h,k2], (103) 

where the coefficient g'gg is given in Eq. (jl92|) . The second amplitude (Fig. [8)3) is a pure massive 
contribution 

M^^^igg ^ x) = Y,iCo{ml,mg)Ti[T''T'^]c^;^^e[fi,i^,h,k2], q = {u,d} (104) 
9 
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Study of the branching ratios of the axi-Higgs. We analyze the dependence on the free parameters 



with u = {u, c, t} and d = {d, s, b}, and the coefficients c^'"? are defined in relations ([88]). The decay 
rate is then given by 

2 
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(105) 



while the expression of the isolated contribution from the corresponding WZ counterterm is instead 
given by 

,3 

(106) 



The decay x ~^ 1^ 
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Figure 3: Branching ratios for the CP-odd scalar Aq of the MSSM 



The partial decay rate computed from the corresponding WZ counterterm and fermion loop, anal- 
ogously to Fig. [H is 
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{Z = Z, Z') which is well defined only for a mass of the Z boson under the threshold production 
mz < rriy-. The couplings are defined as 



X./ ZJ Y f 

= 9z9v c ■ 



(108) 



The vector and axial- vector couplings of the Z bosons to the fermions in the physical basis are related 
to the charges of the chiral fermions by the expressions 



Jv = 2^Qz + Qz )' 9a = -j\Qz - Qz ) 



(109) 



which are obtained, as detailed in [15^ I38j. starting from the interaction basis {W^ , Y, B) by means of 
the following rotations 

9zQ¥ = gvY'^'^O^z + gsY^'fO^z 

gzQ"-/ = 92T''"^0^^2 + gYY'''f0^z+gBY^'^0^z. (110) 

The yj'/-^, Y^/^ and T^^ are the generators of the gauge group of the model in the chiral basis. 
The pseudoscalar triangle Co(?t7-^, m^, mj) involved in the decay x ~^ 7^ with both external lines on 
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their mass-shell, kf = and = m^, is given by (see [50] ) 

Co{ml,m%,mf) = — [m'^^Co{m'^^,mf) - m%Co{m%,mf)] , (111) 

where the structure of Co(m^,m/) has already been studied in ([93} [95]) . In complete analogy, the 
function CQ{w?2,,mf) can be obtained from CQ{m^,mf) just by replacing the first argument with 
m^- Then, the study of the decay rate is closely related to the behaviour of the three-point function 
piip in the various physical domains of its definition. In the domain < mz < my. < 2m f the 
expressions for Co{m'^,mf) and Co(m^,mj) can be read from eq. ([55]) . in particular we obtain 

Co{m%,mf) = ^„ arctan^ — , (112) 

with 



As grows we can have two possible cases. If < mz < 2m f < m^, while the function Co(m^, rraj) 
develops real and imaginary part as shown in eq. (j95|) . the function Co(m^,"i/) is still well defined. 
But finally if < 2m j < mz < m^ also Co{m'^,mf) develops real and imaginary parts, in particular 

Co(m|,m/) = ReCo(m|,m/) + iImCo(m|, m/), (114) 

in analogy to eq. ([95]) . The massless WZ contribution to the decay rate is 

We just remark that in the calculation of T[x ~^ 77) and T{x ~^ 1^) we have neglected the contribu- 
tions coming from the loops generated by the scalar Hq and h^. 

In Fig. [21 for a given value of the Stueckelberg mass Mi = 1 TeV, we study the dependence on 
the free parameter tan/3 = {10,40} and on gs = {0.2,0.6}. The dependence on tan/3 strongly affects 
the branching ratio for the decay into a cc pair, which appears to be suppressed for a large value of 
tan /3 (tan (3 = 40). The plots clearly show the presence of 3 different main regions in which the decay 
channels of the axi-Higgs are rather different. In the region < m^ < 2.8 GeV the dominant decay is 
in the ss and /i/U channels, with a sizeable gluon channel which becomes very relevant around m^ = 3 
GeV. 

For 2.8 < m^^ < 8.5 GeV the dominant decay channel is the rf, followed by a third region with 
> 8.5 GeV in which the hb channel opens up. The 4 plots describe different charge assignments. 
One can notice rather straightforwardly that the leading behaviour in each mass region remains the 
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Figure 4: Study of the leptonic and the quarks branching ratios of the axi-Higgs. We analyze the dependence 
on the function , q^^,Aq^). 



same in each plot, while the subleading channels get reshuffled in their separate contributions. We 
show in Fig. [3] for a comparison, the branching ratios for the CP-odd scalar of the MSSM as a function 
of its mass. In this case the dominant regions are two, divided approximately into the two regions by 
= 5 GeV and the where the dominant decays are into ss (in the lower region) and into rf (in the 
higher mass region). 

• Total rates and dependence on the charge assignments 

We show in Fig. (|4]) plots which illustrate the behaviour of the (inclusive) branching ratios of the 
axi-higgs into quarks and leptons as a function of the mass of the physical axion, obtained by varying 
the charge assignments of the model. The enhancement of the lepton decay channels for a light axion 
mass between 4 and 8 GeV, respect to the quark channel, is very stable against these variations. 
These changes are described by the function f{qQ^,q!^^,Aq^) and in the various cases are are almost 
coincident, and this is due to the fact that the differences in the smaller than 10~^. 

7 CP-even sector: decays and associated production 

We now move to discuss the CP-even sector of the model which involves the two states Hq and /iq. 
We include all the relevant channels, such as the //, the WW, ZZ and Z7 and diphoton channels. 

• Decays into // 



32 



We start by calculating the tree level decay rate into fermions, which is given by 

T{h/H ^ //) = '^\c^lHJ\^N,{f) il - iji^ (116) 

where the scalar couplings to the fermions c^l^'f have been defined in Eqs. (jl20lll9p . The decay for 
the SM scalar Higgs is obtained from Eq. (|116|) just by substituting the coupling c^l^^S with the SM 
one, that is — mj/f, where v is the vacuum expecation value of the SM Higgs field {v ~ 246 GeV). 

• Tree level decays of the scalar Higgs bosons into and ZZ 

The tree level contributions to the total decay rate of the two Higgs due to the decay into a pair 
and a ZZ pair are computed similarly. These are found to be relevant in the case of Hq for a mass 
m//o greater than 100 GeV. In particular we have added the contributions due to Hq/Hq Z*Z and 
Ho/ho W*W that could be significant when the mass of the scalar is close to the thresholds for 
ZZ and WW pair production. 

For the case of a ZZ pair we obtain 



T{H ZZ) = < 



yH _ (f^HO/hO g2 \ (j 40 2 I imsl,\ m.HF{Mz/mH) \f Mr. < m rr < 9 Mr^ 
'-z*z-['-^zz T^Ji;) \'-T^n, + —w) Qom? it Mz <mH < 2,lVlz 



C?T) ' 1^ (3 + ^ - ^) + Tz^z if > 2M^ 



where the coupling C^^^^^ has been defined in the previous sections and — AA4^'^ / tti'jj . Sw:^w 
short notations for sinOw , cos 9w respectively. 
For the case of two charged VF's we have 

2 



T{H WW) = < 



^w*w = {c^T'l^) M^F{Mw/mH) if Mw<mH< 2Mw 

) ' ^ (3 + 4 - ^) + r^.H. if > 2M^ 



Here the coefficient n is equal to 3 if W* — > tb is not allowed, while is equal to 4 if W* tb is allowed. 
Again, we have defined the coefficient = 4M^/m^. 

In the region 1/2 < a; < 1 the function F{x) is defined as follows 

/47 13 1 \ 

F{x) = -|1 - — - — + ^ + 3(1 - 6x2 + 4^4^| ^j^^^^ji 



(117) 



3(1 - 8x2 _^ 20x^) 1 /3x2 - 1 
H — ^ , ^cos"-^ ' 



V4x2 - 1 V 2x3 

Two photon decay of the scalar Higgs bosons h^,H^ — > 77 
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The computation of the decay rate of a CP-even scalar of the MLSOM into a pair of photons is similar 
to that of the SM. It includes the contribution of the spin 1/2 particles (the fermion loop), of the spin 
1 (the W loop) and the spin {H^ loop) and it is given by 
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E ^c{f)Qf-—{-2)rf [1 + (1 - y)f{y)] 
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WW 
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w - 

2 7^2 



(118) 



where H represents Hq or h^, tj = Arrij j m,jj , t^, = 4M^^ / and the function / (r) has been defined 
previously. The scalar couplings of the lighter Higgs boson to the fermions are shown in >CYuk(^) 
and their expressions are 
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V smp Vd V cosp 

while the scalar couplings of the heavier Higgs boson to the fermions are shown in £.Yuk{H) and 
are given by 



rriu mu cos a 
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(120) 



Vu V sin/3 Vd V cos/3 

Here we have used the relations for the expectation values Vu = vsin/3 and Vd = vcos/3 to express 
these couplings in terms of the couplings of the Higgs boson o the SM. The calculation of the rate into 
gluons is similar but we have only the fermion loop 
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T{H^gg) = 



5127r3 
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• Z7 decay of the scalar Higgs bosons 

The last contribution that we consider in the computation of the total decay rate of Hq / ho is the 
decay into Z^. Also in this case we include only the contribution of the fermion loop and of the spin-1 
loop and we neglect the contribution coming from other loops of scalars 
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4mj/M| and Xw = 4M^/M|, while the functions /i^2 are given in [5T] [52]. We report 



where A/ 
them here for completeness 
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The function /(r) has been defined in a previous section while g{T) is given by 

^Jt — 1 arcsin(l/ y/r) if r > 1 



(123) 
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if r < 1 



It is important to observe that in the first line of Eq. (1122^ we have neglected the contribution to the 
the fermion-boson couplings due to the presence of an extra anomalous U{\). As a matter of fact, in 
our hypothesis (Mi = 1 TeV and qb = 0.1 — 0.2), this contribution is found to be very small and for 
this kind of study these couplings can be considered substantially coincident with those of the SM. 
Finally, the total decay rate for Ho /ho will be given as follows 



tot 



gg 



(124) 



7.1 Numerical results 

We shown in Figs. ([SMI) a comparative study of the branching ratios of the scalars Ho and ho in the 
CP-even sector of the MLSOM and those of the Higgs of the SM. While the Ho and the SM Higgs 
appear to be dominated in their decays by the bb channel only below the WW region, the preferential 
decay of the ho is entirely into this final state for all mass ranges. Both the Ho and the ho appear to 
have a more sizeable decay into rr compared to the SM Higgs. The branching ratio for the decay into 
77 appears to be rather small for the /iq in all the mass range, while the Ho and the SM Higgs show, 
for this channel, a similar behaviour. The two-gluons channel also appears to be more significant 
for both states of the MLSOM compared to the SM Higgs, over the entire mass range, while the cc 
channel appears to be rather suppressed in the case of the ho compared to the SM Higgs. Smaller cc 
rates are also found for the Hq respect to the ordinary Higgs. 



7.2 Associated Production of the CP-even states with vector bosons 

Another possible way of detecting the Higgs at hadron colliders is through its associated production 
with a vector boson. Here we calculate the LO cross section for the Ho/ho associated production with 
a W and Z at the LHC and Tevatron and we have made a comparison with the corresponding rates 
for the ordinary SM Higgs. 
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Figure 5: Study of the branching ratios of the CP-even sector. 
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Figure 6: (Study of the branching ratios of the SM-higgs. 



The partonic cross section can be written as 
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where V represents W or Z, the couplings to the fermions are defined as g^j^ = 2Tj^, Qy = 2Tj^—4Q fs'^ 
for the Z, while 9"^ = gy = for the W. The phase space coefficient is defined as X{x,y,z) = 
{1 — x/z — y/z)"^ — Axy/z^. The total cross section as a function of the mass of the Higgs is given by 
the convolution of the partonic cross section with the PDFs luminosity of the quark-antiquark pair 
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qqbar ^ V + H + X, LHC Vs = 14 TeV qqbar -> V + H + X. Tevatron Vs = 1.96 TeV 
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Figure 7: qq ^ H + V + X at LO a.t the Tevatron. 



produced in the initial state which is given by 

/I dx [ - 1 

X] — fHM^lJ'F,m)fH^{'r/x,iiF,m) + {Hi^ H2} (126) 

q,q 

where /^i?, are the factorization and renormahzation scales and represents the quark probability 
relative to the hadron i^i, etc. We have performed the PDF evolution with CANDIA [53] and we 
have used the set MRST 2001 as input distributions, evolved up to /^i? = = Q. The total cross 
section is given by 

CrLo{'mH,IJ-F,IJ-R) = ^qq{T,HF,fJ'R)^{TS)dT (127) 

J TO 

where tq = {My + mH)'^ / S and S is the center of mass energy of the two incoming hadrons. In Fig. d?]) 
we have shown the plots of the total cross section for the LHC and the Tevatron. In the T^-channel 
the cross section of the SM Higgs is smaller that the similar one of the MLSOM due to the Hq, while 
the same cross section for the ho is more suppressed. A similar behaviour is found both at the LHC 
and at the Tevatron. The cross section in the case of the Z follows a similar pattern in all the three 
cases. 



8 Axi-Higgs production at hadron colliders 

The study of the production of the axi-Higgs at hadron colliders is particularly interesting, especially 
for the possibility of having sizeable branching ratios of the two Higgs Hq and /iq into final state 
axions. Before we come to this study, we pause for some observations regarding the scalar potential 
of the MLSOM, stressing on the similarities and on the differences respect to the 2-Higgs doublets 
Model (2HDM) of type II, which is sufficiently general to describe most of the scalar extensions which 
can be envisioned for LHC applications, and to the potential of the MSSM (see ref. |54j ) . 
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Naturally, the most problematic feature of the 2HDM is the presence of a large number of free 
parameters that affect the possibility of unique and specific predictions, due to the different scenarios 
that may emerge at future experiments in regard to the scalar sector. The MLSOM potential is also 
affected by the same problem. In the case of the MSSM instead, the presence of supersymmetry allows 
some relations between the masses and the couplings and between the mass of the gauge bosons and 
their interaction parameters, which provide further constrains on the allowed parameter space. In the 
scalar sector, in this case, there are only two free parameters, which can be identified with tan /5 and 
with the mass of one of the two Higgs bosons |52j . As a result of this, for instance, in the MSSM, some 
Higgs-to-Higgs decays (see ref. [55]) which are possible in the MLSOM, are avoided. Other features 
of the CP-odd sector of the MLSOM are, for instance, the independence of the mass of the axi-Higgs 
from the parameters of the CP-even sector and the existence of a sum rule relating Hq and ho with 
the vector bosons (V), which is also typical of the 2HDM 



8.1 Axion-like interactions 

As we have discussed above, in the MLSOM the specific feature of the CP-odd sector is the presence of 
axion-like interactions which are not found in the 2IIDM and which are the true novelty of the entire 
construction. It is important to remark that while in models containing CP-odd scalars effective 
interactions such as ^o77 induced by the fermion loops are indeed present, they turn out to be 
proportional to the mass of the fermion running in the loop. This mass-dependence, obviously, is 
completely absent in the MLSOM, since the origin of the Wess Zumino terms, which provide these 
interactions, is related to the restoration of the gauge symmetry of the anomalous effective theory and 
not to a mechanism of symmetry breaking. 

In complete analogy to the case of the SM Higgs, the most relevant sector to look for in the 
production of an axion-like particle is the gluon-gluon fusion channel. It is important to point out 
that given the presence of free parameters that are involved in the generation of its mass appearing 
in the PQ-breaking potential, the axion can be searched for in different kinematical domains because 
the model allows both a very light axion with a mass of the order of 1 GeV or less, and a heavier one. 
As stated before, the particular features of the scalar potential render the predictions of the MLSOM 
different respect to the general 2HDM, due to the presence of the b field, and this of course imposes 
some differences in the treatment of the experimental constraints on the allowed parameter space. 

8.2 The parameters 

The free parameters of the scalar potential can be identified by the coefficients {\uu, ^dd, ^ud, ^'ud) 
are contained in the PQ potential and by {bi, Ai, A2, A3), that are contained in the PQ-breaking poten- 
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Figure 8: The two contributions to the gg ^ x production channel. 

tial. The other free parameters are the ratio of the higgs vevs, identified with tan (3, the Stueckelberg 
mass Ml and the coupling constant gs- 

We start our analysis by considering a scenario in which the mass of the Z boson is exactly re- 
produced at Mz = 91.1876 GcV and the bounds on the mass of the extra Z' arc required to be 
compatible with the current Tevatron data. These conditions can be obtained by fixing the value of 
the anomalous coupling gB ~ 0.1, the value of Vu ~ 246 GeV, the value of the Stueckelberg mass Mi 
in the TeV range and tan/3 = 40. These requirements induce also a small mixing parameter between 
Z and Z' (below 10~^), which is also in agreement with current data. Thus, the mass of the particles 
of the scalar sector are identified by the eight parameters listed above. The value of the mass of the 
axi-Higgs is completely governed by the PQ-breaking sector of the potential and one can always find 
a combination of its parameters so that the axion is very light. The other parameters enter in the 
structure of the mass of the two neutral Higgs and the eigenvalues are found to be very sensitive to 
the selection of these parameters. In our case, these have been chosen as follows: 
{\i,X2,h,bi,Xuu,Xdd,Kd} = {-9 10-5,-1 10-6,-1 10-5,5 10-3,6 10-2,5,0.9}, and we have ob- 
tained the following values for the masses of the CP-even and the CP-odd sectors: 
{mHo ~ 122,m^fQ « 15, « 5} (GeV). 

8.3 The invariant mass distribution 

To quantify the cross section of the processes that we are considering, we introduce the invariant 
mass distributions that must be convoluted with the gluon luminosity in order to obtain predictions 
at hadron level. In general, the total cross section for each process can be determined by using the 
following factorization formula 

a{S,^il,^il)= f dii ^\2gi^l,^il)gi^2,|J^lMc^si^iR),Q^/^^R,Q^/l^F) (129) 

Jo Jo 

where r = Q'^/S, 5(^2 ^/^j^') is the gluon density, function of the Bjorkcn variable ^ and of the factor- 
ization scale fi£,\ A similar expression holds for the invariant mass distributions for the production of 
a pseudoscalar with an invariant mass Q, which is given at parton level by 
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Here |A^p represents the square of the matrix element for the production of n scalar particles in 
the final state, the variables 1^1,^2 represent the fraction of the momentum carried by the partons 
in the collision and d^n is the Lorentz invariant phase space. The invariant mass is defined as 
s + t + u = Q^, while the fraction l/Q^ is the partonic flux. Then we can write at hadron level 

where 

E (132) 

pol,spin 

and the gluon luminosity is given by the following convolution product 

^gg{x,^ll) = f %iy,^l)g(^,f,l). (133) 

Jx y y 

The computation of this cross section for the production of the axi-Higgs pp ^ gg ^ x ~^ ^ via 
gluon fusion involves two contributions: the fermion loop correction and the direct (contact) decay 
due to the Wess-Zumino term, as shown in Fig. [8l with the WZ counterterm suppressed as 1/Mi and 
therefore quite subleading respect to the first. 

At parton level the production cross section for the axi-Higgs via gluon fusion is related to the 
decay rate by the following relation 

(^gg^xi^) = ^\r2 ^(x ^ 99)S{s - m^) = a°gg^J{s - m^) (134) 

where s is the squared partonic cm. energy and A'^c = 8 is the color factor for the gluons. At hadron 
level the total cross section for the inclusive axi-Higgs production is given by 



1 

a{pp -^gg^X + X)= / dT^gg{T)agg^^{TS) = -^(^gg^x^gg(.T)\r. 

J mils ^ 



Q2 



imys 

(135) 

where the variables S and stand for the squared cm. energy of the incoming hadrons and the 

invariant mass of the gluon pair, respectively. In Figs. [9] we show the plots of the total cross section 
at LO at the LHC and at the Tevatron respectively, for the production of the axion and of each of the 
CP-even Hq, ho Higgs, and the corresponding plots for the SM Higgs. Notice that the result shows a 
ratehr sharp rise of the production cross section with a decrease of the axion mass, larger by a factor 
of 10 compared to the case of other CP-even scalars. A similar rise is found also for the CP-odd sector 
of the 2HDM, being typical of the CP-odd sector. 
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Figure 9: Cross section for tlie production of the two Higgs ho and Hq and the axi-Higgs via gluon-gluon fusion 
at LO at the LHC (left panel) and at the Tevatron (right panel). 

9 Axion plus photon production 

In this section we compute the production of an axion plus one photon at the LHC in leading order 
(LO), given by the diagrams in FigllOl The computation of the amplitude requires the three-point 
correlator between two photon and one axion, with one off-shell photon and with mj ^ 0. This can 
be achieved by using the parametrization of the trilinear vertex with two off-shell external legs and 
away from the chiral limit (see |26j). 

Denoting by T^'^'^ the correlator with outgoing momenta , k!^ and incoming momentum k^ , the 
generalized WI gives the following relation 



where the tensor T^"^ is defined by 



iCo{k^,klm^.)m} 

-e[ki,k2,n,u\ 



(136) 



(137) 



Performing the change of the momenta ki ki,k2 — > —q,k k2, we obtain the expression for 
the three-point correlator between two photons (one off-shell) and one pseudoscalar, suitable for our 
calculation. The function Co has the following expression 



where we have defined 

a2 = 



02 + 1 
02 — 1 



log^ 



2/03 + 1 



03 



1 



03 



Arnj 



(138) 



(139) 



We can identify four kinematic regions in which the function Cq can be analytically continued: 
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Figure 10: Production channel for a single axion plus a photon 
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The squared and averaged partonic contributions are given by 



spin 



rt2 + ^^2 



m^2 ,mf,, . 
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spin 
1X /,,.. ;^ /^X 



(140) 
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mr, 

Qfe2g^^Y.Q}e^Re[Cois,m^2,m},)] (144) 



where 0^,/vf' is Ofi/fu for an up type quark, while 02i/vd for a down type quark. 
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Figure 11: Invariant mass distribution for the associated production of an axion plus one photon at the LHC. 



Integrating over the two-particle phase space we obtain 
1 1 (s - m2)2(s + ml) 



CTi(s,m^) 



487riV2 2i 



Qje 



/' 



'^Q'j,e'^—j-Co{s,m^2,mj,) 



2 ^"^/' 



awz[s,m 



1 1 (s-m^)^(g + m^) 
72 '^/^ (^77) ' 



487riV2 2s 
1 



487rA^2 2s 



1 (s - m^)2(s + m^) , Of, mi 



(145) 



where dint denotes the interference term. Introducing the invariant mass distribution at hadron level, 
we have 



da _ o-(Q^m^; 
dQ2 ~ s 



(146) 



where the parton luminosity has been previously defined and Q represents the invariant mass of 
the final state. 

We show in Fig. [11] a plot of the cross section for the production of an axion and one photon at 
the LHC as a function of the mass of the x- The mass dependence of the result is quite small, except 
for a larger mass of the particle, in a region where it is Higgs-like. For an ultralight axion the value 
of the cross section is around 10"^ pb. We have shown the contribution from the triangle and the 
Wess-Zumino terms combined and separately, in order to show the dominance of one channel respect 
to the other. The Wess-Zumino term is indeed strongly suppressed (by a factor of 10^^). 
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Figure 12: Channels for multi axion production from gluon gluon fusion. 

10 Multi axion production 

One of the peculiarities of a light axion-like particle is its possibility to generate cascade decays with 
multi- lepton final states which are more sizeable especially for a mass of x ™- the GeV range. We 
have indeed seen that for around few GeV's, the largest contribution to the branching ratio of 
its decay is predominantly into leptons, and for this reason we are going to investigate systematically 
this particular interval in parameter sace. Our analysis will include two types of vertices, the trilinear 
XxHq, ho vertex and the vertex. As we are going to see, multilepton decays will be sizeable even in 
the presence of a considerable phase space suppression and we will quantify them rather accurately. 

We consider both the production of axions in combination with a scalar of the CP-even sector of the 
MLSOM, and final states made entirely of several light axions which branch primarily into leptons. We 
consider the gluon fusion channel, in which the production of the CP-even scalars (/iq, Hq) is mediated 
by the top and bottom loops. The sizeable values of the multi-axion cross sections for the invariant 
mass distributions are related to the large production cross sections which are typical of pseudoscalar 
channels and to the large values of the reduced couplings - normalized to the SM ones - of the trilinear 
interactions of the scalars. The leading contribution to the production cross section comes from the 
fermion loop graph with a final state axion. In the model, each contribution is accompanied by the 
corresponding WZ counterterm, which is suppressed by a factor of 10^ compared to the loop graph 
(see Fig. E]). 

Channels involving several final state axions can be built rather easily, list of several diagrams 
contributing to these channels is given in Fig. [T2j For instance, the simplest process involves a gg — ho 
production channel combined with the ho — XX vertex. In this case the WZ counterterm is absent. A 



44 



similar process is the gg — x triangle vertex, followed by the xx^o vertex, which gives the combination 
of a X and of a CP-even Higgs [ho) in the final state. In this case the channel is accompanied by 
a WZ term gg — x describing the direct interaction of the two gluons of the initial state with the 
axion. Cascading channels can be easily obtained by combining trilinear (xx^o) and quadrilinear (x^) 
vertices, which are more involved and that we will study below. 

10.1 i/o, ho — 3x decay 

The amplitude for the on-shell production of three x and one scalar Higgs - through the process 
gg — > hPxXX " is given by the sum of a part containing the fermion triangle plus the counterterm 

M = Mloop + Mcount ■ (147) 
Defining s = = {ki + k2)^ we can write the square of the matrix element as 

2^ \Mh,3x\ - \ (JV2_1)2 2^(4 ) C^^hi99sCx^:;^^c\f{rf)\ + 8/^2 _ 1)2 

spin,pol. f \ c J 



(iV2-l)2 ttV j (g2-m|)2(l-a;2 + ^)2' 

(148) 

where X2 = p' ■ q/q^ ■ The coefficient r/ is defined as Anrj/q'^, while pi = 4m^/g^, and p2 = Arrv^/q^. 
The coefficient g^GG of the counterterm is defined in the previous sections and the couplings of the 
axion to the up-type and down-type quarks are given by 



c. 



uu _ X _ "^uVd dd_ rn.^_m^ 



The details of the computation can be found in an appendix. 
10.2 4 - X decay 

We move to discuss the possibility of producing four axions in the final state mediated by a CP-even 
higgs (i^o, ^o)- At parton level, the squared amplitude for the process gg ^ B. ^ Arx^s, given by 

2. 1-^4x1 - (^2 _ 1)2 2. (^2 _ ^2 )2(i _ ^^)2 |1 + (1 " r/)/(r/)| (150) 

spin,pol. f \ 

where H = Hq, ho and the couplings of the higgs to the quarks are given by 

uu _ '^U jj uu _ p 

Ch — -0-12 Cfj — U22 

V V 

i?i2 = — cosct i?22=sinQ;. (151) 
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Figure 13: gg 2-scalar reactions mediated by trilinear vertices 



The coefficients Ri2,R22 are the matrix elements for passing from the interaction eigenstate basis to 
the physical basis, already defined in the previous sections. 

The plots for the production of four scalar particles via gluon-gluon fusion are shown in Fig. [TH 
Notice that the production of four axions and that of three axions and one ho show invariant mass 
distribution which are rather similar in their sizes. This is due to the fact that in this study we have 
chosen ho to be not too much heavier than x ("^/lo ~ 15 GeV). Details on the computation of the 
4-particle phase space can be found in an appendix. We have performed a direct computations of 
the phase space integrals, which have been reduced into a 2-dimensional form and then have been 
integrated numerically. The results of this study are shown in Fig. [13] for the Tevatron and the LHC 
respectively. The plots presented in the two figures show sizeable rates which become large on the 
Higgs (Hq) resonance, chosen to be at 120 GeV. At the LHC the peak value of the cross section for 
PP ~^ XXi niediated by the Hq is larger by a factor of about 10 — 100 compared to the Tevatron and 
would be significant. In the same figures, the same production channel, mediated by the /iq, is also 
resonant at 15 GeV, but is not shown in our study since it involves an extrapolation of the parton 
distributions towards the small-x region, which we have not included in our analysis. 

Coming to the 4-axions final state, the numerical values of the various distributions are shown in 
Fig. [m where they appear to be down by a factor of approximately 10^ compared to the analogous 
ones with 2 x's or with one x aiid one CP-even Higgs in the final state. We have summarized in 
Tab.dSD the numerical value of the cross sections at a representative value of Q at which they appear 
to be sizeable, within the parametric choices used in our analysis. The largest values shown are those 
on the resonances of the two neutral Higgs. The multilepton channels, for a GeV axion, appear to be 
rather small even on the largest production resonance, which is on the peak of the Hq, due to a large 
phase space suppression. Typical resonant rates are 10~^ pb/GeV for 4 muons and 10~^^ pb/GeV 
for the production of 8 muons. For final states with 8 muons mediated by the ho in the non-resonant 
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Figure 14: Production of 3 and 4 scalars from gluon gluon fusion mediated by trilinear and quadrilinear 
vertices. 



region and coming from the pairwise decays of 4 axions, the rates are much smaUer (~ 10 GeV 
pb/GeV). 



Process 


Q 


da/dQ (LHC) 


Q 


da/dQ (T) 


gg ^ ho^ 4x 


45 


^ 10-'^ 


22 


4 • 10-4 


gg ^ Ho^ 4x 


Mho 


103 


Mho 


1.56 


99 ^ X^^X + ho 


50 


5 • 10-4 


40 


2 • 10-^ 


99 ^ X ^ ^X + Ho 


150 


2 • lO"'^ 


150 


f^lO-8 


gg ^ ho^ 2x 


45 


26 


20 


2.5 • 10^ 


gg ^ Ho^ 2x 


Mho 


324 • 10^ 


20 


4.9 • 10^ 


99 ^ X ^ ho + X 


45 


0.69 


20 


2.5 • 10^ 


gg ^ x^ Ho + x 


150 


« 10-3 


150 


« 10-5 


gg ^ Ho ho + ho ^ 4x 




5 • 10^ 


150 


82 



Table 8: A list of processes analyzed at hadron colliders at the LHC and at the Tevatron (T). Q is in GeV and 
da/dQ in pb/GeV. 



11 The light mass region of the axion and its hfetime 

One obvious question to ask is whether the axi-Higgs, which takes the role of a valid example of a 
gauged axion, has any chance of being a dark matter candidate, with properties which remain quite 
distinct from those of the axion of the PQ model. 

As we have already remarked in the introduction, axion-like particles originate from the gauging 
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Axion lifetime, gB= 0.1, M,= lTeV, f(-l,-l,4) 
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Figure 15: Lifetime for an ultralight axion as a function of tan/3 



of anomalous symmetries, and take the role of phases in the scalar potential, which is characterized 
by a small curvature in these variables. We have seen that for a particle mass in the GeV range the 
branching ratios for its decay into leptons appear to be too large for the particle to be long lived. We 
can pause for a moment and try to understand the origin of this result. 

The axion interaction with the fermions is generated by the ordinary Yukawa couplings, being the 
particle part of the scalar sector of the model. In particular, the CP-odd contributions are re-expressed 
in the physical basis by the elements of the rotation matrix together with other parameters, the 
most significant of them being /3, as shown by Eq. (j85p . Notice that the matrix elements of this matrix 
are 0(1), which means that we can't expect a large suppression of its coupling to the fermions just 
from its mixing with the other CP-odd components of the Higgs sector. 

If we look more closely into the two contributions which appear in the decay of an axion, the 
triangle diagram and the WZ term, one finds that the contribution from the triangle is 0{mf/v), 
where v is the vev which represents the symmetry breaking scale of the symmetry to which the axion 
is associated as a phase of a complex scalar. Consider, for instance, the mechanism of chiral decoupling, 
that we have described in the previous sections. In this case, the only interaction of the axion with 
the gauge fields takes place through the WZ terms, since there are no Yukawa couplings between the 
light fermions and the pseudoscalar. Then, if we assume that the decoupling scale M^i = gsMg is 
around lO^'^ GeV, which is the decoupling scale of a right handed neutrino in a typical leptogenesis 
scenario, the decay rate is simply given by the relation 




(152) 
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where g^^y and ggg are proportional to 

« ^ (153) 

and is dominated by the 2-photons and 2-gluons channels. For a very weakly coupled axion, with a 
small value of the coupling constant {qb « 10~^), we have indeed a long lived particle of around 1 
GeV with a rather long lifetime 

r^ = ^^ 10^6 s. (154) 

In the MLSOM instead, the suppression comes from the Stiickelberg mass Mi while the Yukawa 
couplings remain unsuppressed. Therefore, in this model, the structure of the axion-fermion-fermion 
interaction is proportional to mj /v x , where v is of the order of the electroweak scale and is 
of order 1 if Mi is in the TeV region. In these conditions, the MLSOM allows a long lived axion only 
if this is very light, with a mass rriy^ « 10~^ eV, which is again, specific of this construction. 

We show in Fig. [15] plots of the hfetime of a very light axion (10~^ — 10^^ eV) of the MLSOM as 
a function of tan /3, which shows that in both cases the particle is very long lived, with features which 
resemble quite closely those of the traditional Peccei-Quinn axion. 



12 Conclusions 

In this work we have an analyzed the phenomenology of the physical axion that emerges in several 
extensions of the Standard Model and which include an anomalous U{1) gauge symmetry. We have 
focused our study on a mass window characterized by an axion of a light-to- intermediate mass, which 
is probably easier to detect at colliders, although windows for a particle of even lower mass can be 
analyzed in a similar fashion. One of the most appealing features of the class of models that we have 
presented consists in the possibility to justify in a natural way a particle in the CP-odd sector of such 
a small mass, which would be more difficult to motivate at theoretical level in other constructions. We 
have shown that the origins of the class of effective actions that are characterized by the presence of 
such a state could be quite different. For instance, in the case of brane models, the small mass of the 
axions is parameterized by extra terms in the potential which are identified by the symmetry of the 
low energy model and in which the axion appears as a complex phase. These terms may induce a small 
tilting on the scalar potential, giving a small mass to the physical axion, extracted after electroweak 
symmetry breaking. A similar tilting is induced by the instanton vacuum in the case of the Peccei- 
Quinn axion, and as such, it is possible, given the strong analogy between our case and the PQ case, 
to borrow most of the results - well known in the case of the invisible axion model - and extend them 
to this more general model. A very light axion would be, with no doubt, a good candidate for dark 
matter. 
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We have also shown, although in a simplified model, that effective actions which resemble quite 
closely the MLSOM, can be obtained by a completely different approach, using the decoupling of a 
chiral fermion - due to a large vev of a Higgs to which this fermion is coupled - from the effective 
theory. The charge assignments of generalizations of the MLSOM can be obtained by this approach. 

In this second case our analysis has to be considered rather preliminar and needs further extensions, 
although we expect that most of the features of the special form of chiral decoupling that wc have 
proposed can be worked out more closely in the context of a Grand Unified Theory. The generalization 
of this analysis to the supcrsymmctric case appears to be rather interesting as arc the cosmological 
implications of the presence of a gauged axion (with or without supersymmetry) in the low energy 
spectra of these theories which deserve further studies. 
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13 Appendix A: The Lagrangean 

The classical lagrangean of the model is explicitly given by 

Co= - [F^^F^n - IttI F^^F^n - \f^,F^^'^ - \fI,F^>^^ 

+ \{d,. + i92'-^Wi+igyq'^Al + igj,'^B^)Ha\^ 

+ QLi ir {d^ + ig^^Gl + ig^'-^Wl + i5y(??"^< + igBQ^'^B^ Qu 

+ urn i^i" (d^ + igyq^'^'^Al + igsq^^'^^B^) um + dm i^i" (d^ + igyq^"^ Al + igBQ^B^^B^) dm 

+ Li ir {d^ + ig2^W/, + igyq^^Al + igsQ^^^B^ L, 

+ ejii i-fi' (d^ + igyq^'^'^Al + igBq^^^B^ em + il" (d^ + «5y9?''^^^ + igBQB^^B^ vm 

+ ]^{d^h + MstB^f 

+ V{H^,Ha,b), (155) 

which generates Sq. We have summed over SU{3) index a = 1, 2, 8, over the SU{2) index j = 1,2, 
and over the fermion index i = 1,2,3 denoting a given generation. We have denoted with tl 



the 
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field-strength for the gluons and with the field strength of the weak gauge bosons TV^. and 
F^j, are the field-strengths related to the abelian hypercharge and the extra abelian gauge boson, B, 
which has anomalous interactions with a typical generation of the Standard Model. The fermions are 
either left-handed or right-handed Dirac spinors fi, /r and they fall in the usual SU{3)c and SU{2)\y 
representations of the Standard Model. 



14 Appendix B. Matrices of the potential 



The mass matrix in the CP-even sector is given by 

A/'2(l,2) = 2 (Sv^Ag cos^ /3 + 3v^A2 sin^ /3 + 2u2Ai sin2^ - 2i;2A„<isin2/3 + 6) 
A/'2(2,2) = -2sec/3(-4Add?;^cos^^-3A3i;^sin/3cos^/3-|-A2?^^sin^/3-|-6sin/3) 



-2(-4v^ Au„ sin^ P + v'^Xs cos^ /? cot /? - -^v'^M sin 2/3-1-6 cot /?) 



In the CP-odd sector we have 



/ cot/3 



A/3 = --VuVdC^' 



1 



-1 

tan P 



Vd 



VuVd 



Mi 

ll-i'd 
Ml 



In the charged sector, the mass matrix elements are 

A/Kl,!) = -2cot/3(A3Cos2/3+ (Ai - A'„d)sin2/3 + A2sin2/3)i;2 -26cot^ 

A/"i(l,2) = 2 (A3 cos^ /? + (Ai - A'„d) sin 2/3 A2 sin^ /?) + 26 

7Vi(2,2) = -2(A3Cos^/3-F(Ai-A'„d)sin2/3 + A2sin^/3)f2tan/3-26tan/3. 



(156) 



(157) 



(158) 



15 Appendix C. Matrix and quadrilinear interactions 

We report for completeness the matrix O^, which is given by 
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^ — N = -NcosP (159) 



" VuVd 
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^1 '^\IW^^^d 
1 



Vd 
° I'm I'd 



—N = Nsm.p (160) 



51 



31 



12 



22 



32 



Ml 

-(nB-nB 



= NQi cos /3 



sin/3 



-.2 



U 2 + 2 



cos/3 



(161) 

(162) 

(163) 
(164) 



(Ox; 



13 



(Ox; 



23 



33 



1 + 



Ml 



AT 



(gf - ) 

Ml 
1 



Vu COS /3 



1 + 



(gg-gf )^ vivl 



COS /3 = iVQi cos (3 

{<lu-1d)\ 



-N 



Ml 



u„ cos /3 
= = N.. 



Ml 



sin/3 = -iV(5isin/3 



The coefficients appearing in tfie quadrilinear vertices are given by 



Ri """" = sin a cos a [iddiO^if - luu{0 
Rt""""" = \ cos' a [iddiO^if + luuiO^if] 



[iddiOlf + luuiO'^if] 

= sin a cos [{Olf - {0^,f] 



R^ 



RxW^_^^^2^l^^ [(0|i)2 + (0 



j^x'H^o ^ sinacosa/i [{O^if - {O^if] - Asina cos ah^^O^i (Oji + O^i) + 0{1/M^) 



Aq 



B 



(165) 



(166) 
(167) 



Ml 



2^0^0 1 
Iti — — — 

3 2 



cos^ aZi [(0^i)2 + {O^if + AO^iO^ 
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21 J 



+ cos2 ah^ [vd{^0\,Ol, + 202^ Of,) - v^{AO^,Ol^ + 20^,03^)] + 0(1/Mf) 
= -\ sin^ «Zi [(O^^)^ + {Ol,f - ^0\,0\\ 

+ sin^ ah^ bd(-40^iOfi + 20fi03\) - v^{-AO^,Ol^ + 20^iOfi)] + 0{l/Ml 

A -B 

Rf"°''° = sin a cos aO^^O^^ik - h) + sin a cos a— ^Of^ [vdOr^ih - 3/3) + ■y«02i(^3 - 3/2)] 

iVL\ 

RfH^H^ = i cos^ a {/2 [O^iO?, - (0?,)2] + /3 [O^iO^i - {0^^?] } 

+ cos^ a^Ofi {^^u [0^1 («3 + 3^2) + 2/20^J - Vd [0^1(^2 + 3/3) + 2/302^]} 



(168) 



+ sin^ a^Ofi [0^1 (l3 + 3^2) - 2^20^1] - Vd + - 2/302^]} 



Appendix C: Axi-Higgs Trilinear Interactions 

RfH' = cos a [{O^^fvdldd - {O't^fvuluu] , 
Rf"' = sin a [{O^ifvdldd + {O^ifvJuu] , 
Rt"° = cosalud [{O^ifvu - [Otifvd] , 
Rf''" =-smalud m^fv^ + {0^,?^ , 

Rf'^°=WsinaOl^{0^,-0^,), 
Rf"" = cosah [0^,(20^, + 0^,)vu - 0^,(01 + 202^)^^] 

+2cosa/i^Ofi [O^^Vdivd - 2vu) + O^^Vu(.Vu - 2vd)] + 0{1/M^), 
Rf^° = -small [Ofi(0?i - 202^,)!;^ - 02^(20^1 + Ol)vu] 

= ^ cos a [0^1/3 i20^iVd + O^ivu - Vd)) - O^.h {0^i{vd - Vu) + 20^,Vu)] 
Aa^ 

+ cos a^Ofi [-vuh (30^1 + O^^K - 2vd)) - Vdh {'iO\^Vd + 0^^{vd - 2vu))] , 
Rf^° = ^ smap^^h {20^iVu + OfiK + ^d)) + O^^h {O^iivd + v^) + 20^^Vd)] 

53 



Aa^ 

+2smah-^Ol [O^M'^Vd + v^) - O^M^d + 2vu)] + 0{1/mI), 



2Mi 



31 



(169) 



17 Appendix D: Quadrilinear self interactions in the CP-even sector 



For Hq we have 



while for /iq we have 



2 

For the interactions of the type h^^HQ^ wc obtain 



jj4 . 1 

i?i ° = (cos a) - {Xuu + Add) , 

= -(cosa)^^A„d, 
Rf° = (cosq;)^^Ai, 

i?f° = (cosq)'^^(A2 + A3), (170) 



/j4 1 

Ri" = -(sin a) (A„„ + Add) , 
^2° = -^(sinQ;)^A„d, 

= ^(sina)'^Ai, 

i?J° = J(sina)4(A2 + A3). (171) 



Ri° ° = -(sin a)^ (cos a)^ (A„„ + Add) , 
j^Hihl ^ _(sina)2(cosa)2Ai. (172) 



For the interactions of the type ho Hq we obtain 



Rf°^° = (sina)^ cos a (Add - Ku) , 



while for h^IIfp we obtain 



i^f o'^o = (sin af cos a (A3 - A2) , (173) 



j^oH^ = (cos af sin a (Add - Ku) , 

^hoH^ = _(cos af sin a (A2 - A3) , (174) 
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18 Appendix E: Trilinear self interactions in the CP-even sector 



For we have 



Jj3 „ 

Ri ° = (cos a) {vdXdd - VuXuu), 
R2 ° = (cos ay{vd - Vu)Xud, 
i?3 ° = {cosaYivd - Vu)Xi, 

R? = \{cosaf [i;„(3A2 + A3) - Vd{\2 + BAs)] , (175) 



while for we have 



For the case H'^Hq we have 



i?/ = {sm.ay{yd\dd + VuKu), 
^2° = -{sinayivd + Vu)Xud, 
rI° = {smaf{vd + Vu)Xi, 

nf = ^(sina)3 [z;„(3A2 + A3) + VdiX2 + SAa)] • (176) 



j^hoHo ^ ^cosa{sma)'^{vdXdd - VuXuu), 
= C0SQ;(sinQ:)^Aud('U(i - Vy), 

Rf"° = -\ cosa{s\naf{vu + Vd){X2 - A3), (177) 



while for the case hoHg we have 

j^f^oHo ^ 2 gin q:(cos a)'^{vdXdd + ■"uAuu), 
1^2°^° = sinQ;(cos a)'^Xud{vd + Vu), 

= ^ sina(cosa)2(i;d - Vu){X2 - A3). (178) 

19 Appendix F. The axion Lagrangian in the physical basis 

We have seen that after symmetry breaking, in the scalar sector we isolate a physical axion, Xi 
also called the axi-Higgs. Here we present the axion Lagrangian rotated on the basis of the mass 
eigenstates. In particular, the VF3, and B gauge bosons become linear combinations of the physical 
states A^, Z, Z' . Indeed, the mass-matrix in the neutral gauge sector is given by 

jCmass = {Ws, Y, B) M2 
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where B is the Stiickelberg field and the mass matrix is defined as 



^ . 92^ v'^ -92 9y^'^ -92^3 
M^ = -^-g2gyV^ gy\^ gyX^ | (179) 

-92 9y^b 2Mf + iV^B ^ 



with 

,B2„,2 I „B2„.2\ „2 _ („B„.2 , R„.2 



Nbb = {luK + irvj) 9l Xb = (gf < + q^vj) 9b- (180) 
Here Vu and Vd denote the vevs of the two Higgs fields Hu, Hd while and are the Higgs charges 



under the extra anomalous U{1)b- We have also defined v = yv^ + v"^ and g = \J 92+ 9y- "^^^ 
mass-squared eigenstates of the mass matrix corresponding to one zero mass eigenvalue for the photon 
and two non-zero mass eigenvalues for the Z and for the Z' vector bosons, are respectively given 

by 



^ {2MI + g\^ + Nbb - yJ{2Ml-g-^v^ + Nj,j,f + ^g^xl^ (181) 



5V J_5^, J_5^.^ 2 2^ 



^ 2 

The mass of the Z gauge boson gets corrections of the order /Mi converging to the SM value as 
Ml 00, while the mass of the Z' gauge boson can grow large with Mi. The physical gauge fields 
can be obtained from the rotation matrix 



(183) 



which can be approximated at the first order as 

f + 0(6f) + 0(6?) §61 . (184) 

V -fei ^61 1 + 0(6?)/ 

Moreover, after symmetry breaking, as we have already shown in eq. (jlOp . the Stiickelberg field b is 
rotated by means of the matrix as follows 

b = Olx + 0^,GUO^,Gl, (185) 
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where the elements of the rotation matrix have the following expressions 

0^1 = . ^ 03^3 = ; ^ 0^2 = 0- (186) 

{lS~Q^? vWd ^ V ^ 



Then, starting from eq. (|185p . the Goldstone modes and G^' in the 7-basis are obtained by the 
combination 

Gl = C'zG^ + C'z,G^' . (187) 

More details can be found in [23]. Starting from the WZ Lagrangian in the Y-basis 

C'^^'iZsis = DbTr[F^ AF^]+FbTr[F^ AF^] 

+GYYhF^ A + CsBbF^ A F^ + GyshF^ A F^ (188) 

and rotating into the physical mass eigestates using eqs. (jl83p and (jl85p we obtain the axion-like 
terms of the WZ Lagrangian 

+ g^zz XF^AF^+ x F^' A F^' + g^^ xF^AF^ 

+ g^z'XF'' AF^' +glz,xF^ AF^\ (189) 

^axion^QZ^ = (^gG^Tr[F^ AF^] + cl_G^Tr[F^+ AF^-]+c?;^G^ F-^ AF-^ 
+ 4z F^ AF^ + c%z' F^' A F^' + c^^ G^ F^ A F^ 
+ c^z'G^ F"/ AF^' + 4z'G^ F^ AF^', (190) 

where Z stays for Z, Z' . Finally the WZ Lagrangian in the physical basis is given by the sum of three 
contributions 

where we have identified the physical couplings of the axi-Higgs x to the gauge bosons as 
9l- = FOl 

g^^ = {FO^,,Oi.,^, + CyyO^^O^^) 

9zz = {FO{^^zOw^z + GyyOyzOyz + GbbO^zObz + GybOyzObz) 
9z'Z' = {FOw^z'OwsZ' + GyyOyz'Oyz' + GbbObz'Obz' + GybOyz'Obz') 
g^z = {2FO^^^^O^^z + ^GyyO^^O^z + GybO^^O^z) Ol 
g^z' = {2FO^^^O^.^z' + '^GyyO^^O^z' + GybO^,O^Z') 
9zz' = i^FOwazOwaZ' + 2CyyOyzOyz' + 2GbbObzObz' 

+ GybO^zO^z' + GybO^z'O^z) (192) 
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and the interactions of the NG bosons {Z = Z^Z') with the gauge bosons 



cf z = {^Oy^^zOwzZ + CyyOyzOyz + CbbO^zObz + CybOyzObz) O^^C'z 
^%'z> = i^O^^z'OwaZ' + C'yyOyz'Oyz' + CbbO^z'Obz' + C'ybOy^'O^^/) O^^C'z 
c^z = i^FO^s-yOwsZ + 2CyyO^^O^^ + CybO^^O^z) O^^C'^ 
c^z' = {'^FO^s'y<sZ' + 2CyyO^^O^z' + CysO^^O^^,) Og^gC^ 

cf Z' = (^-^O^^gZ^^TygZ' + SCyyOy^Oyz' + 2Cbb0^zObz' 

+ CybOyzObz' + CybOyz'Obz) O^^C'^- 

(193) 

where Z stays for Z, Z' . We also summarize for convenience the coefficients of the WZ counterterms 

Ml 2 2 ^ ' 
Ml 2 ^ ' 

CbB = J^-^an^BBB^ 
'^YY - J^Wy-^J->BYY^ 

Cyb = ^WY^Dy^B, (194) 
with a„ = — and the chiral asymmetries have been defined, for brevity, in the following way 



/ 



^B — 


f 


^B — 


Q 


F>BBB = 




F>BYY = 


f 


F>YBB = 





Q 



f 



58 



20 Appendix G. Three- and Four-particle phase space 



The three and four body phase space in the case of massive particles can be computed directly in four 
dimensions since there are no soft and collinear divergences. The reactions that we are considering 
are 5(^1) +5(^2) ^ xiff) ^ x{p)x{r)x{p') H{p') and g{ki) + g{k2) H{q^) xip) x{r) x{p') x{p'), 
where the on-shell conditions are given by = r'^ = = m^2 and p'^ = m\ for the first reaction, 
while for the second we have = r'^ = p^ = p'"^ = m^2. The computation follows closely [SH], with 
some modifications due to our specific case, given the three axions and one higgs boson in the final 
state. 

20.1 Phase space for the three axions and one scalar higgs final state 

In four dimensions we can write the most general formula as follows 

d^4=— / , / — -:T—n — r^ri^TrfS^iq-p-p'-r-r'), 196) 

2\ J 2po(27r)3 2p'o(27r)3 2ro(27r)^ 2r^(27r)^^ y w J, \ J 

where 1/2! is a statistical factor that takes into account the fact that a pair of identical particles are 
produced in the final state. The reference frame in the CM of r, r' can be chosen as 

r = (ro, \r\ sin sin 0, \ f\ sin cos 0, \r\ cos 9) 

r = (ro, — |r| sin0sin(/>, — |r| sin0cos(/>, — \f \ cos 9) 

/ TTi^ 

V Po 



P =Po(l>0, Wl -:rsm.a,^l -^rcosa). (197) 

V Po V Po 

We introduce the following variables 

xi = 2^-, X2 = 2^-, y = 2 2 ' 

qz qz 

ro = V^^, 1^ = Jrl-m? 



2 ' 

9 2 
mi, 

pi = 4^, P2 = 4^. 198 
qz qz 

From the momentum conservation equations {q — p)"^ = {p' + r + r') and {q — p')^ = (p' + r + r') we 
derive the expression of po and Pq as a function of the variables (xi, X2, y, \/q^, ni^, mu) as follows 

(l-X2-y)v^ , rn]j-ml 
Po = o ^ — + 



2vV\/y 
2\/y 2v^^ 



p'o = ^'-^'-y^^ + I^ipi. (199) 
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Using the equation = (p + p' + r + r')^, we obtain the expression of cos a in terms of the kinematic 
variables defined above 

Po + P'o + \f(?yjy - {pI - ml) - (p'l - mjj) - 
cosa = = . (200) 

2^pl-ml^p'l-mjj 

In order to integrate the expression given in Eq. (jl96p it is useful to introduce the following identities 
"^'^ i27r)^5\t -r-r') = l, q^ [ -^(27r)5(t2 - q^y) = 1 (201) 



/ 



(27r)4^ ' ' 7 (27r) 

which allow us to incorporate r and r' in the t state. Thus, we obtain 

i / W / 2^ / lii^P-)**'' — ') 



W^^.WTT^^i'^^fnt -r-r'), (202) 



(27r) '""'y 2ro(27r)3 2r^(27r)3 

where &{y) is the Heaviside step function. In this way we have factorized the expression of phase 
space as a product of d^3 x d^2 

^•^3= / TTr4^ / Tr4^ / -^^i2nr6iq-t-p-p'). (203) 



2to(27r)3 7 2po{2TTf J 2^0(2^)3 



Integrating over (i<I>2 we obtain 



2n 



(i$2 = 77^9^/1-— / '^^ / # (204) 



4(27r)2V y 70 

where we have defined v = 1/2(1 — cos0), while the integration over brings us to 

{2tt)_ \p\^d\p\n3 \p' \'^d\p'\ sin P d(3n2 
2to 2po(2vr)3 2p'o(27r) 

where Iq and (3 have been computed below 



Finally we obtain 



*0 = \/\t\'^ + q'^y =\/W + \pT + '^\P\\P'\ COS p + q^y, 

cos/? = ii^--^--2?-M-i-l-P2)-i4-Pi) ^206) 
2Y/rcf - P2y/x'^ - pi 

d^3 = J dxidx2, (207) 

and the final result for the d^4 phase space is given by 

d$4 = — ^— / l-L±dy / dxi / dx2 dv #, (208) 



2!(4vr)6y^2 V y -^-2- 

where the integration limits are discussed in the next section. 
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20.2 Integration limits 

At this stage we need to define the integration hmits of the integrals appearing in the four body phase 
space. From the definitions of xi,X2 it is clear that < xi,a;2 < 1, but imposing the reahty condition 
of the square root we obtain 

xi>^/p^ X2 > VpT- (209) 

Solving the condition — 1 < cos/? < 1 with respect to X2 we obtain a bound on this variable which is 
given by 

= _ ^ 4^ {[(^1 - 2)(pi - 4(xi + y - 1)) + (xi - 2)p2] 

±\/ (xl - P2)[l6xj + 8x1 (4y + /?i - /O2 - 4) + I6y^ + (pa - pi + 4)^ - 8y{pi + p2 + 4)]| . 

(210) 

Again, we have to impose the condition X2+ > ^/pl which gives us a condition on the variable xi 

My) < — ^ r— , 211) 

4-2^ 

but xi must be such that the square root in Eq. (|210p is real 

xi{y) < \{-^y - 4^/pr^/y - /oi + /02 + 4) 

XI > ^2. (212) 
From these three conditions we can extrapolate some conditions on the y variable 

Pi<y<\{^x^ ^2-1? (213) 

20.3 Phase space for a four axions final state 

In the case of a four axions final state we have a simplification in the computation since P2= P\ = P- 
Thus, the four body phase space is computed exactly as in [56j and the final result is given by 

where the factor 1/4! is a statistical factor that takes into account the four identical particles in the 
final state and the integration bounds are defined as 

1 

X2± 



(2 - xi)(2 + p-2y- 2xi) ± 2 J(x2 - p)[{xi -l + yf- py] 



4(l-xi) + p 
xi < 1 - y - ^/py 

yi+ = (1 - ^pf. (215) 
These integrals have been computed numerically. 
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